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The  problem  of  a  mixture  of  two  stimulants  in  a  biological  quantal 
assay  is  investigated  from  a  mathematical  standpoint.  The  basic  assump¬ 
tion  is  made  that  the  response  region  does  not  depend  on  biological  con¬ 
siderations  -  i.e.,  given  a  specified  mixture  of  stimulants  z  ,  the  response 
region  is  defined  by  the  point  z'  in  the  p-variate  space  where  there  are 
p  stimulants  under  consideration;  instead,  the  probability  functions, 
themselves,  may  take  on  different  forms.  A  general  form  is  proposed  and 
investigated.  Three  analytic  models  (one  utilizing  the  bivariate  normal 
distribution,  one  a  bivariate  logistic  distribution  developed  by  Gurtbel 
(1961) ,  and  one  a  bivariate  Burr  distribution  developed  by  this  author) 
are  employed  in  this  investigation.  The  investigation  includes  the  anal¬ 
ysis  of  data,  under  the  three  analytic  models,  which  had  been  classified 
by  previous  investigators  as  examples  of  synergistic  action,  simple  similar 
action,  independent  action,  and  additive  action.  The  residual  analyses 
are  included  as  well  as  the  FORTRAN  IV  subroutines  used  in  evaluating  the 
functions,  the  partial  derivatives  and  the  weights. 

The  investigation  lends  some  support  to  the  assumption  of  a  constant 
response  region  for  a  diversity  of  mixtures  of  stimulants.  The  analytic 


iv 


model  incorporating  the  bivariate  Burr  distribution  is  recommended  for 

a  1.x  taaca  uiilcas  Ll ic  of  pilTwITtC^C  “3  tC  tC  £  £>  tli.IT.2  ^ C  ^  1.C  2  p2Ti!T.2.2Ty 

concern,  in  which  case  the  analytic  model  utilizing  the  bivariate  normal 
distribution  is  recommended.  The  bivariate  Burr  distribution  developed 
ir.  this  paper  is  found  to  be  more  useful  in  application  than  that  devel¬ 
oped  by  Takahasi  (1965)  . 
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CHAPTER  I 


1.  Introduction 

The  joint  action  of  mixtures  of  stimulants  in  a  biological  assay 
has  been  investigated  by  Bliss  (1939),  Finney  (1942),  Plackett  and  Hewlett 
(1967),  Ashford  and  smith  (1966),  and  others.  Plackett  and  Hewlett  have 
made  their  investigations  largely  from  the  standpoint  of  biological  con¬ 
siderations  such  as  the  p.iysiology  of  the  biological  organism  being  used 
in  experimentation.  Ashford  and  Smith,  on  the  other  hand,  have  dealt  with 


the  problem  somewhat  more  within  a  mathematical  framework.  In  this  paper, 
the  problem  will  be  approached  mathematically. 

For  the  purposes  of  this  paper  a  biological  assay  of  a  mixture  of 
two  stimulants  will  be  conducted  as  follows:  A  population  of  N  organisms 
is  divided  at  random  into  t  groups,  where  the  i  group  is  of  size  n^  , 
n^  +  n2  +  • •  •  +  nfc  =  N  .  The  i^1  group  receives  a  treatment  of  a  pre¬ 
determined  mixture  (z  ,  z  )  of  two  drugs,  where  z,.  is  the  quantity 

11  *1  31 

of  stimulant  j  measured  in  any  convenient  units,  r^  is  the  observed  number 
which  manifest  a  prescribed  quantal  response.  The  observed  relative 
frequency  of  response  p^  =  r^Ah  is  an  estimate  of  the  probability  of  an 
organism  responding  if  picked  at  random  from  the  population.  The  proba¬ 
bility  that  this  organism  picked  at  random  will  respond  when  treated  by 


the  mixture  (z^  ,  may  be  assumed  to  take  on  a  general  form,  say, 

P(Zli  '  Z2i  '  • 

Now  the  probability  of  r^  responses  with  the  ifc^  combination  of 


1 


levels  of  drugs  can  be  written  as 


f  Kt^i 


r 


cLTCn,-xLU  Lr'"li  '  *2i  '  Lx  "  rv"li  '  ‘2i  '  *'J 


..iV'i 


r  -  0,  1,  2,  •  ••  ,  t 


elsewhere . 


A  series  of  t  combinations  of  doses  is  tested  in  an  experiment.  The  prob¬ 
ability  of  a  particular  set  of  r^'s  is  equal  to  e:tp(L)  ,  the  likelihood, 
whore 


t  t  t 

L  «  l  r  In  (P.)  +  l  (ni~ri)ln(2i)  +  £  lnln^l/r^l  (n^r^  l]  (2) 

i«l  i=l  i=l 


and  P.  =  P(z,,  ,  z_,  ,  0) ,  Q,  =  1  -  P.  .  The  maximum  likelihood  estimator 

i  11  21  -  l  1 

0  of  a  parameter  0  ,  0  an  element  of  0  ,  must  satisfy  the  relation 


t  r.  3P.  t  (n.-r.)  3Q,  t  n.(p.-P.)  9P. 

„  3L  r  i  _ i  .  p  i  i  *i  r  i  i  i  _ i 

'  »8  ^  »t  30  h  Si  30  Ji  PA  30  • 


Direct  solution  for  0  is  not  in  general  possible,  but  iterative  techniques 
are  available  which  give  a  convergent  series  of  approximations  to  the 
solution. 

The  following  procedure  for  two  parameters  0  and  ,  0  and  (f>  elements 

of  0  ,  is  of  completely  general  applicability  and  may  easily  be  extended 

for  the  estimation  of  a  greater  number  of  parameters.  By  the  Taylor- 

3  L  9L 

Maclaurin  expansion  of  —  ,  —  (See  relation  (3).)  ignoring  quantities 

30  09 

containing  terms  of  higher  than  the  first  degree 


3L_  sJl  .  52L 

301  d0  302  %  39^ 


(4) 


3 


34> 


3^  0 


=  0 


where  the  addition  of  the  suffix  .1  to  0  f  4>  indicated  that  the  first 
approximations  are  to  be  substituted  after  differentiation.  The  solutions 

6  ,  6  are  adjustments  to  0,  ,  i .  which  give  the  improved  approximations 

0  4>  ±  i 


°2  ■  8I  +  *0  ’  *2 


’1  +  5* 


Equations  (4)  may  be  simplified  through  the  following  procedure 
which  will  be  illustrated  by  means  of  the  first  of  equations  (4). 


or 


_3L  +  ,  afi+  6  .sV- . 

301  50  a02  %  80!**! 


.  3^L  32L  =  _3L 

0  8Q2  "  ♦  9013*1  “  901 


0  30 


W,  _3_  f  3L~]  3L 

%  3*!  1^3©!  J  "  3®! 


-  6 


0  30i 


-  6 


<j>  3  <t>. 


nipi 

!!i. 

t 

I 

n, 

X 

3pi 

P ,  Q. 
1*1 

30! 

Oi 

30x 

n.p, 
i  i 

3P._ 

t 

I 

i=l 

n. 

i 

9fi 

P.Q, 

li 

* — 1 
© 
CO 

«i 

3©1 

-  6, 


(  t  n.Pi  32P. 

v 

ffiV  s  Vi 

^3pi  V2 

U-ipA  ^ 

i“i 

i-i  plq; 

^901/ 

t  n,  32P.  t  n.  /3P. \2 ) 
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(  t  n.p.  3  P.  t  n.p.  3P,  3P 

)  y  _ i_  _  y  x1  .  i  _ 

l  <■  v  n  an  a*  c  ■>  ao  a* 


2, 


c 

niPi  <)t'i 

i  _  y 

“i  a  r. 

i 

i-1 

p  o?  30i  3* 

1  L 

L 

1  i=l 

Ql  3013d>1 

t 

y 

ni  3Pi  api  i 

l  mJk 

L 

i=l 

o?  30i  a*i : 

f  301 

At  this  stage  the  equation  may  be  simplified  by  putting  p.  =  P  in 

32L  32t 

the  coefficients  of  — -  ,  ~  %  ,  i.e.#  on  the  left  hand  side  of 

30*  3013*1 

the  last  equation  to  give  expected  instead  of  empirical  values.  The  last 
equation  then  reduces  to 


The  latter  equation  in  (4)  can  be  reduced  by  means 
i.e.,  putting  =  P^  in  the  coefficients  of 
Thus  equations  (4)  are  simplified  to 


of  a  similar  procedure, 


32L  32L 

3 *2  '  ^ 


1 


t  n.jp^,^  /3P.\ 

iil  Pil2il  \30l/ 


(5) 


t  n.lp.-P.^ 


I 

i=l 


PilCil 


Here  the  addition  of  the  suffix  1  to  P..  ,  Q..  indicates  that  the  first 

il  ll 

approximations  are  to  be  used  in  the  evaluation  of  P  (z^  ,  z^  ,0)  . 
Equations  (5)  illustrate  that  only  first  derivatives  are  needed  in  this 


iterative  procedure. 
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Now, 


u.. 

•"l 


2.  Methodology  for  Obtaining  Estimates 
it  will  be  seen  from  ’bat  follows  that  relation  (3)  can  be 
linear  least  squares  [Moore  and  Zeigler 


(1967)].  Assume  that  the  data  corresponds  to  the  mathematical  model 


y^  *  ,  o)  +  ,  i  =  1,  2,  •••  ,  t 

where  the  y.  are  observed  random  variables,  is  a  vector  of  known 
independent  variables,  a  is  a  vector  of  unknown  parameters,  and  ci  is  a 
random  variable  such  that  E(e^)  =  0  ,  E(e^)  =  c^  ,  and  E(e^e_.)  -  0  for 
all  i  ft  j  .  Then  the  vector  of  unknown  parameters  may  be  estimated  by 
minimizing  the  weighted  sum  of  squares , 


s  =  l  (yL  -  h(zi  ,  a))\  ,  <7) 

i=l  ' 

where  W.  is  an  appropriate  weight.  If  the  usual  procedure  is  modified  so 
that  the  partial  derivatives  are  taken  ignoring  Wh  ,  the  normal  equations 
are 

t  3h(z.  ,  a) 

-21.  =  -  2  T  W.  [y.  -  h(z.  ,  a)]  - r- - -  0  ,  (8) 

3\  i=l  1  1  _1  '  ^ 

£or  ^  =  1.,  2,  ,  S.  ,  where  l  is  the  number  of  unknown  parameters.  Now 

by  letting  W±  =  ni/pi2i  (the  reciprocal  of  the  variance  of  pi)  ,  Y±  =  Pi  * 
h(z.  ,  a)  =  P.  ,  Q.  =  1  -  ,  and  =  ©  it  can  be  seen  that  relation 

(3)  and  equation  (0)  are  equivalent.  Thus  the  maximum  likelihood  estimate 
can  be  obtained  by  means  of  a  modified  weighted  non-linear  least  squares. 
Relations  (5)  and  their  equivalent  extensions  are  used  in  the  modified 
non-linear  least  squares  fitting  of  equation  (6) . 
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3.  Consideration  of  Necessary  Conditions  on  P(z^  ,  z2^  ' 
- r—1 - i - - 

From  this  point  onward  the  vector  z.  =  (  I  will  he  considered 

v 

from  the  standpoint  of  a  mixture  of  stimulants  where  a  transformation  has 
been  applied  to  the  original  dosage  levels  so  that  -»  is  equivalent  to 
zero  dosage  and  +°°  is  equivalent  to  an  infinite  dosage.  For  the  purposes 
of  this  paper  P  (z^  ,  z^  ,  0)  must  satisfy  the  following  conditions 


P<Zli  ' 

,  0)  =  P0«  ,  z2i  ,  0)  =  1  , 

(9) 

'<v  >  ~ 

/  0)  “  P1(Zli  t  0^)  9 

(10) 

and 

P<~  ,  z2i 

,  0)  =  P2(z2i  ,  02)  , 

(ID 

where  P^Cz.^  >  0.)  and  p2^Z2i  '  §2^  are  not  *n  9eneral  zero,  but  rather 
are  marginal  probabilities,  i.e.,  the  probability  of  a  random  individual 
biological  organism  responding  if  it  is  given  a  dosage  of  stimulant  j 
corresponding  to  .  Conditions  (9) ,  (10) ,  and  (11)  imply  the  conditions 


p(_oo  ,  -0°  ,  0)  =  0  , 

(12) 

and 

P(+“  ,  +»«  ,  0)  =  1 

(13) 

All  five  of  these  conditions  are  necessary  in  a  bioassay  of  quantal  response 
data  involving  a  mixture  of  two  stimulants.  Natural  extensions  of  these 
conditions  for  a  mixture  of  more  than  two  stimulants  are  now  obvious. 
Plackett  and  Hewlett  (1967)  proposed  that 


f(zu  ,  z2.)dz1.dz2i 


(14) 


where  f(z^  «  z^)  is  a  bivariate  density  with  the  usual  properties  and 


R  is  defined  on  the  basis  of  biological  considerations,  thus  implying  that 
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the  region  of  integration  may,  arbitrarily,  be  changed  due  to  biological 
considerations.  Their  papers  do  not  indicate  any  homogeneity  in  the  regions. 
Nowhere  is  there  a  aeneral  formulation  for  P(z_.  ,  z. .  ,  G)  where  the  form 

J.J. 

of  the  region  of  integration  is  homogeneous,  much  less  constant.  It  would 
appear  that  the  region  of  integration  should  be  constant  except  possibly 
for  simple  monotonic  transformations  of  the  original  dosage  levels,  such 
as  a  logarithmic  transformation.  The  bivariate  function  itself  might  be, 
in  specific  instances,  of  different  types  but  still  retaining  a  constant 
response  region. 

Let  f’1(z1  ,  0^) ,  F2  (z2  ,  G25  be  univariate  distributions  where  the 
parameter  vectors  0^  ,  ©2  are  not,  in  general,  equal.  Note  that  F^z^  >  t 

F2*z2  '  ®2^  are  not  necessarily  even  from  the  same  family  of  distributions, 
e.g.,  the  family  of  normal  distributions.  Let  F,(z1  ,  z2  ,  0)  be  a  bivariate 
distribution  such  that  F^z^  >  +“  /  0)  ”  Fi^zi  '  and  F3(+“  »  z2  ,  §2) 

=  F2(z2  ,  ©2 ) ,  where  O’  =  (0^  ,  §2  ,  Gj)  .  Now,  what  is  needed  is  a 
function  which  satisfies  conditions  (9)  through  (13) . 

Let 

H(z1  ,  z2  ,  0)  =  F1<z1  ,  0  )  +  F2(z2  ,  02)  -  (z^  ,  z2  ,  0)  .  (15) 

Then  H(-«  ,  z2  ,  0)  =  F2(z2  ,  02> ,  H(z1  ,  -»  ,  0)  -  F]L(z1  ,  ©x) , 

H  (+°°  ,  z2  ,  0)  =  1  =  H(z1  ,  +®  ,  0),  H(-°>  ,  ,  0)  =  0  ,  and  H(-H”  ,  +*  ,  0) 

=  1  .  Thus  H(z.^  ,  z2  ,  0)  does  satisfy  conditions  (  9)  through  (13)  which 
suggests  that 

P(«1±  »  z2i  •  0)  =  Fi(zi  '  +  F2(z2  '  -2}  "  F3(Z1  '  z2  '  (16> 

is  a  general  formulation  for  Pfz^  #  Z2i  '  where  the  response  region 
is  constant.  Note  that  the  forms  F^z  ^  ,  0^  ,  F2  (z2i  ,  02>  ,  F^z^  ,  z2i  ,  0) 


8 

are  completely  general  distributions  whose  forms  can  depend  on  biological, 
i  .  _  _j_i _ _ _ &+■  4-v»o  cam*  +* i mt>  H\f»  recrion  of  inte ~ 

uieiiuwaif  01  wwubj.  .  *‘-  — •  -  — 

gration  is  constant  and  easily  understood  from  a  geometrical  standpoint 
as  well  as  from  other  standpoints.  It  is  noted  here  that  the  general 
formulation  for  P  (z  ,  z21  ,  0)  can  easily  be  extended  for  a  vector  of 
more  than  two  stimulants.  The  utility  of  this  form  is  quite  general? 
the  only  restrictions  being  conditions  (9)  through  (13)  which  have  been 
imposed  in  the  development  of  the  general  form  in  equation  (16) . 


CHAPTER  II 


It  is  natural  in  the  study  of  a  mixture  of  two  stimulants  to  con¬ 
sider  a  bivariate  probit  or  normit.  Probit  analysis  has  no  advantage  over 
a  normit  analysis  if  the  analysis  is  run  on  a  high-speed  computer.  Also, 
the  analyses  are  equivalent.  Almost  all  of  the  work  that  has  been  done 
to  date  has  been  along  the  lines  of  a  bivariate  normit. 

Bliss  (1939)  was  among  the  first  to  study  the  action  of  mixtures 
of  two  stimulants.  He  classified  the  joint  action  of  two  stimulants  into 
three  biological  categories!  independent  joint  action,  similar  joint 
action,  and  synergistic  action.  Independent  joint  action  occurs  when¬ 
ever  two  components  act  on  different  vital  systems  in  the  organism  and 
do  not  interact  with  one  another.  Similar  joint  action  is  observed  when¬ 
ever  two  components  act  independently  of  one  another  but  on  the  same  vital 
system.  Synergistic  action  is  characterized  by  a  larger  frequency  of 
response  than  could  be  predicted  from  experiments  using  the  individual 
stimulants.  He  mentions  antagonistic  action  but  did  not  treat  this  concept 
at  all.  He  stated  that  it  is  the  reverse  of  synergistic  action.  Finney 
(1942)  suggested  that  antagonism  is  negative  synergism  and  can  thus  be 
treated  in  the  same  category  as  synergism. 

Bliss  (1939)  ,  for  the  category  of  independent  joint  action,  plotted 
expected  response  in  probits  against  dosage  of  mixture  in  logarithms.  At 
each  point  the  ratio  of  the  amount  of  a  given  stimulant  to  the  amount  of 
the  other  stimulant  was  held  constant.  These  curves  were  not  smooth  but 
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rather  fell  into  two  segments  each  of  which  appeared  to  be  a  straight 
line.  The  transition  from  one  straight  line  to  the  other  was  relatively 
abrupt.  He  suggested  the  equation 

Pc  =  PA  +  PBd  "  PA>  U  “  r)  (17) 

where  pA  >  PB  f  PA  i-s  the  probability  of  response  due  to  the  effect  of 
stimulant  A  ,  p  is  the  probability  of  response  due  to  stimulant  B  ,  r 
is  a  measure  of  "association  of  susceptibilities,"  and  pc  is  the  proba¬ 
bility  of  death  due  to  the  combination  of  stimulants  A  and  B  He  did  not 
indicate  what,  if  any,  relation  he  assumed  between  equation  (17)  and  the 
plots  of  data. 

For  the  category  of  similar  joint  action  Bliss  suggested  equation 

y  =  a'  +  b  log(D  +  kD.)  (18) 

C  Ad 

for  the  dosage  response  curve,  where  and  Dfi  are  the  respective  doses 
of  stimulants  A  and  B  in  the  mixture  and  k  is  the  ratio  of  the  frequency 
of  response  of  the  individual  stimulants.  The  plots  for  this  case  are 
thus  straight  lines. 

Bliss  suggested  two  possible  equations  for  synergistic  action.  The 
first,  which  relates  the  total  amount  of  active  material  (Dft  +  DQ)  and  the 
amount  of  the  more  active  stimulant,  say  A  ,  is 

<Da  +  DB)Dj  =  k  ,  (19) 

where  D  and  D  are  in  original  dosage  units,  which  implies  the  probability 
A  B 

of  response  to  the  combination  of  the  two  stimulants  is  determined  by  the 
sum  of  the  ingredients  multiplied  by  some  power  of  the  amount  of  the  more 
active  stimulant.  The  second  equation,  again  with  A  being  the  more  active 


stimulant,  is 


11 


(1  +  k  D  )dJ  =  fc  (20) 

1  A  B  2 

which  was  suggested  for  the  cases  where  the  proportion  of  A  approaches 
zero.  It  should  be  noted  that  these  suggested  equations  do  not  bear 
any  clear  logical  relation  one  to  another. 

Plackett  and  Hewlett  (1961)  utilize  the  following  biological 
classification  of  joint  drug  actions: 


Similar 

Dissimilar 

Non-interactive 

Simple  Similar 

Independent 

Interactive 

Complex  Similar 

Dependent 

Here  the  suggestion  is  that  the  actions  of  the  stimulants  are  similar  or 
dissimilar  respectively  as  the  stimulants  act  on  the  same  biological  site 
or  on  different  ones,  and  as  interactive  or  non-interactive  depending  on 
the  presence  or  absence  of  synergism  (or  antagonism).  They,  then,  propose 
mathematical  equations  (some  in  an  implicit  form)  based  on  the  above 
biological  classifications  which,  again,  do  not  bear  any  clear  logical 
relation  one  to  another.  Finally,  they  introduce  a  statistical  concept 
into  their  presentation  by  making  an  assumption  as  to  the  bivariate 
distribution  of  ,  z where  z^  ,  z are  the  respective  tolerances  to 
stimulants  A  and  B  .  They  suggested  that  a  reasonable  assumption  would 
be  that  the  log  tolerances  log  z^  ,  log  z are  distributed  bivariate 
normally.  They  did  not  give  examples  of  data  fitting  any  of  the  proposed 
models . 

Ashford  and  Smith  (1964)  approached  the  problem  somewhat  differently. 
They  classified  the  mathematical  model  as  interactive  or  non-interactive 
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ratlier  than  attempting  to  classify  on  the  basis  of  biological  considera¬ 
tions.  Thev  define  non-interaction  as  being  equivalent  to  the  condition 
on  P  =  P(z^  >  x2  '  *  the  Probability  of  response,  where  z^  and  z are 

the  logarithms  of  dose,  such  that 


2  2 

W  (P)  =  PP(PP  -PP  )  +  p  (P  P  -  p  p  \  =  0 
12 1  J  12'  2  12  1  122;  12'  1  22  2  ll' 


3P  32P  93P 

where  P  =  t —  ,  P  0  ~  r — 5 —  ,  and  P  „  =  r — - — 5 —  .  Their  mathematical 

a  3z  afi  3z  3z„  ctBy  3z  3z„3z 

a  a  (3  a  g  y 

classification  is  not  equivalent  to  Plackett  and  Hewlett's.  Ashford  and 


Smith  remarked  that  no  valid  distinction  can  be  made  between  similar  and 


dissimilar  action  purely  on  the  basis  of  quantal  response  data. 

Ashford  and  Smith  published  some  trivariate  data  on  exposure  to 
coal  dust  for  which  the  response  was  the  prevalence  of  pneumoconiosis  for 
groups  of  mine  workers.  The  three  dosage  variables,  respectively,  were 
the  time  spent  in  years  at  coalface  coal-getting,  coalface  preparation, 
and  elsewhere  underground.  They  assumed  that  the  tolerances  were  normally 
distributed.  They  then  compared  two  models  where  the  regions  of  response 
were  not  only  different  but  were  each  complicated  functions  of  the  dosage 
levels.  They  applied  chi-square  goodness-of-fit  tests  (each  with  fifteen 
degrees  of  freedom)  to  the  models  obtaining  chi-square  values  of  12.73  and 
16.86  ,  respectively,  from  which  they  quote  the  corresponding  approximate 
significance  levels.  They  do  not  indicate  explicitly  the  form  of  the 
probability  function  used  but  rather  only  the  functional  forms  indicating 
the  response  regions. 

Zeigler  and  Moore  (1966)  presented  a  paper  at  the  126+-b  Annual  Meeting 
of  the  American  Statistical  Association  on  "Multivariate  Quantal  Response 
Analysis  Using  Regression  Methods."  In  this  paper,  in  addition  to  showing 
that  weighted  least  squares  can  be  used  to  converge  on  maximum  likelihood 
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estimates,  they  fitted  a  bivaric.be  normal  distribution  to  toxicity  data 
involving  the  direct  sprays  of  Pyrethrins  and  D.D.T.  in  Shell  Oil  P31 
apipxied  to  flour  beetles  (Tribolium  castaneum) .  Using  a  chi-square 
goodness-of-fl t  test  with  nineteen  degrees  of  freedom,  they  obtained  a 
value  of  12.17  and  reached  the  conclusion  that  the  fit  was  satisfactory. 

None  of  the  investigations  up  to  this  point  have  utilized  the  general 
form  suggested  in  Chapter  I,  although  the  specific  form  utilized  by  Zeigler 
and  Moore  (1966)  is  equivalent  for  the  special  case  where  the  tolerances 
z and  to  drugs  A  and  B  are  each  distributed  normally. 

It  would  seem  useful  to  do  some  numerical  studies  utilizing  some 
of  the  data  in  the  literature  with  some  analytic  models  which  conform  to 
the  general  form  in  equation  (16) .  For  this  purpose,  seven  sets  of  data 
were  utilized.  Included  among  these  were  sets  that  have  been  classified 
in  the  following  categories  by  previous  investigators:  synergistic  action, 
simple  similar  action,  independent  action,  and  additive  action. 

Data  set  one,  classified  as  synergistic  by  Bliss  (1939),  was  first 
published  by  Kagy  and  Richardson  (1936) .  This  set  is  from  a  study  of  the 
combined  action  of  2-4-dinitro-6-cyclohexylphenol  and  petroleum  oil 
sprayed  in  emulsions  against  eggs  of  a  plant  bug  (Lygueus  kalmii  Stll . ) . 
Data  set  two,  published  by  Plackett  and  Hewlett  (1952)  ,  was  classified  by 
them  as  simple  similar  action.  This  data  set  is  from  a  study  of  the 
combined  action  of  D.D.T.  and  methoxychlor  applied  in  Shell  Oil  P31  to 
flour  beetles.  Data  set  three,  published  by  Hewlett  and  Plackett  (1950), 
was  classified  by  them  as  independent  action.  This  is  the  data  set  which 
Zeigler  and  Moore  (1966)  fitted  to  a  bivariate  normal  by  means  of  weighted 
least  squares.  Data  sets  four,  five,  and  six,  published  by  Martin  (1942), 
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were  not  classified  by  the  investigator  into  any  category.  Data  set 
four  is  from  »  study  of  the  toxicity  of  the  combined  action  of  rotenone 
and  a  dequelin  concentrate  in  a  medium  of  0.5%  saponin  containing  5%  of 
alcohol  applied  to  chrysanthemum  aphides  (Macrosiphonie  11a  sanbomi) . 

Data  set  five  is  from  a  study  of  the  toxicity  of  the  combined  action  of 
rotenone  and  H-alliptone  under  the  same  laboratory  conditions  as  data  set 
four.  Data  set  six  is  from  a  study  of  the  toxicity  of  the  combined  action 
of  rotenone  and  i-a-toxicarol  under  like  laboratory  conditions.  These 
three  data  sets  showed  some  signs  of  synergism  to  the  investigator,  but 
he  did  not  find  it  to  be  significant  in  any  one  of  the  data  sets.  Data 
set  seven,  published  by  Ashford  and  Smith  (1964) ,  is  from  a  study  of  the 
prevalence  of  pneumoconiosis  in  groups  of  mine  workers  where  the  years 
spent  on  "coal-getting"  is  one  imput  and  the  other  imput  is  years  spent 
in  "haulage."  This  data  set  was  classified  as  an  example  of  additive 
action  by  the  investigators. 

A  bivariate  normit  analysis  was  run  on  the  above  seven  sets  of  data. 
The  analytic  model  for  the  bivariate  normit  analysis  was 

/*al+BlZl  -  — 

,  z2  ,  0)  =  l  (2tt)  2  exp (-  —  t2)dt 


VB2Z2 


_  1 

(2ir )  2  exp(-  —  s2)ds 


(22) 


-VB1Z1  ,VB2Z2 


I 

"  —00  "  —co 


(2tt/T'-  pZ' ) 


.-1 


exp[-(t2  -  2pts  +  s2)/2(1  -  p2)]dtds,  -«  <  z  , 


z_  <  » 
2 


15 


A  modified  least  squares  (see  Chapter  I)  FORTRAN  IV  Computer  program  was 
utilized  on  a  Model  44  IBM  360  system.  A  resume  of  the  results  is  given 
in  u'anie  i. 

The  following  is  a  brief  explanation  of  the  items  listed  in  Table  .1 
as  well  as  the  next  two  tables:  N  is  the  number  of  stimulant  combinations. 
SSE  is  the  weighted  sum  of  squares  due  to  error  which  is  approximately  dis¬ 
tributed  as  a  chi-square.  SSR  is  the  weighted  sum  of  squares  due  to  re¬ 
gression  and  is  computed  as  SST  -  SSE  where  SST  is  the  weighted  sum  of 

squares  adjusted  for  the  weighted  mean.  SSR  is  approximately  distributed 
2 

as  a  chi-square.  R  ,  which  is  computed  as  SSR/SST  ,  tells  what  portion 

2 

of  SST  is  due  to  regression.  Computing  SSR  as  SST  -  SSE  and  R  as  SSR/SST 

gives  both  a  conservative  estimate  of  the  significance  of  regression  and 

2 

a  conservative  coefficient  of  determination  R  .  The  column  entitled 
"No.  of  significant  chi-squares"  tells  how  many  of  the  chi-square  statistics 
computed  at  each  dosage  level  (stimulant  combination)  exceeded  3.84  ,  the 
.95  value  of  a  chi-square  with  one  degree  of  freedom. 


Data 

Set 

No. 

N 

No.  of 
Significant 
Chi-squares 

SSE 

d.f. 

SSR  | 

1 

d.f. 

R2 

1 

18 

5 

67.029 

13 

66041  | 

4 

.99899 

2 

10 

1 

21.775 

5 

598.86  | 

4 

.96491 

3 

24 

0 

11.805 

19 

11176  | 

4 

.99894 

4 

17 

2 

27.147 

12 

1656.0  | 

4 

.98387 

5 

12 

2 

28.947 

7 

921.36  | 

4 

.96954 

6 

15 

0 

10.145 

10 

30672  | 

4 

.99967 

7 

40 

2 

38.141 

35 

217.75  | 

4 

.85095 

TABLE  1 
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For  all  of  the  data  sets  the  regression  is  found  to  be  significant 
using  SSR  as  the  indicator.  However,  the  chi-square  for  departure  from 
the  model  is  insignificant  in  only  three  of  the  cases,  namely  data  sets 
three,  six,  and  seven,  which  include  the  cases  of  independent  action  and 
additive  action. 

The  synergistic  data  (data  set  1)  had  sample  sizes  ranging  from 
240  to  479  (see  Appendix  I)  at  its  eighteen  data  points.  The  bivariate 
nonnit  analysis  indicated  that  five  of  these  points  differed  significantly 
from  the  bivariate  normal  model.  Some  of  these  points  were  marginal  data 
points  and  some  were  not.  One  of  the  data  points  contributed  34.266  to 
the  cumulative  chi-square,  slightly  more  than  half  of  the  total,  but  the 
chi-square  would  still  be  significant  even  without  this  particular  data 
point.  Upon  examination  of  the  residuals,  the  fit  does  look  good  with  the 
exception  of  the  one  data  point,  but  with  the  large  sample  size  at  each 
point,  the  fit  would  have  to  be  extremely  close  in  order  for  the  cumulative 
chi-square  to  be  insignificant.  On  the  whole,  it  is  felt  that  the  bivariate 
nonnit  analysis  did  quite  well  with  the  data  and  that  the  model  does  des¬ 
cribe  the  phenomenon  reasonably  well,  considering  the  significance  of 

regression  (SSR) ,  the  weighted  sum  of  squares  due  to  error  (SSE) ,  along 

2 

with  the  sample  sizes,  and  the  coefficient  of  determination  R 

The  simple  similar  action  data  (data  set  2)  had  sample  .sizes  ranging 
from  148  to  200  (see  Appendix  II)  at  its  ten  data  points.  The  analysis 
indicated  that  one  of  these  points  differed  significantly  from  the  bivariate 
normal  model.  Again  upon  examination  of  the  residuals,  the  fit  does  look 
good  although  not  quite  as  good  as  the  previous  data  set.  The  conclusion 
based  on  the  analysis  of  the  data  is  that  the  model  does  describe  the 
phenomenon  fairly  well,  with  the  -xception  of  the  one  data  point. 
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The  independent  action  data  (data  set  3)  had  sample  sizes  ranging 
from  48  to  50  (see  Appendix  III)  at  its  twenty  four  data  points.  The 
model  does  fit  the  data  well  and  none  of  the  data  points  differed 
significantly  from  the  model.  The  weiqhted  sum  of  squares  due  to  error 
is  11.805  .  Zeigler  and  Moore  (1966)  fitted  this  same  data  set  and  the 
weighted  sum  of  squares  due  to  error  for  their  model  is  12.17  ,  thus 
indicating  the  similarity  of  the  fit. 

Data  sets  four  and  five  are  quite  similar.  They  had  sample  sizes 
ranging  from  28  to  51  (see  Appendices  IV  and  V)  at  their  data  points. 

Each  had  two  data  points  that  differed  significantly  from  the  bivariate 
normal  model  and  examination  of  the  residuals  does  not  indicate  as  good 
a  fit  as  for  any  of  the  previous  data  sets.  The  model  still  does  describe 
most  of  the  data  points  well,  but  it  does  not  seem  to  do  as  well  as  for 
the  earlier  cases. 

Data  set  six  had  sample  sizes  ranging  from  48  to  51  (see  Appendix  VI) 
at  its  fifteen  data  points.  The  model  does  fit  the  data  well  and  none  of 
the  data  points  exhibit  a  significant  deviation  from  the  model.  Two 
bivariate  normit  analyses  were  run  on  this  data  set  using  slightly  different 
convergence  criteria.  The  first  run  utilized  the  relative  change  in  the 
unweighted  sum  of  squares  due  to  error  and  the  second  the  relative  change 
in  the  weighted  sum  of  squares  due  to  error.  The  first  run  after  conver¬ 
gence  had  the  sum  of  squares  due  to  error  as  0.031265,  while  the  weighted 
sum  of  squares  due  to  error  was  0.78159  x  1015  .  The  second  run  after 
convergence  had  the  sum  of  squares  due  to  error  as  0.032040  while  the 
weighted  sum  of  squares  due  to  error  was  10.145  .  Which  criteria  produces 
the  best  fit  becomes  questionable  at  this  point.  It  would  seem  that  either 
set  of  parameter  estimates  would  have  to  be  considered  acceptable  despite 
the  large  chi-square  value  attributed  to  the  first  fit. 
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Data  set  seven  had  sample  sizes  ranging  from  2  to  135  (see  Appendix  VII) 
at  its  forty  data  points.  The  model  does  fit  the  data  well  although  there 
are  two  data  points  which  deviate  significantly  from  the  model.  Ashford 
and  Smith  (1964),  who  classified  this  data  as  an  example  of  additive  ac<  on, 
fitted  the  data  to  a  model,  assuming  the  marginals  to  be  logistic,  using 
a  rather  complicated  response  region  which  does  not  seem  to  have  been 
necessary. 

In  general  the  bivariate  normit  analysis  seems  to  do  quite  well  with 
a  diversity  of  mixtures  of  stimulants,  as  is  evidenced  by  the  seven  sets 
of  data  analyses  here.  These  analyses  appear  to  lend  support  to  the 
assumption  that  the  form  of  the  response  region  should  remain  constant 
irrespective  of  the  biological  consideration.*- ,  at  least  in  relation  to  a 


bivariate  normit. 


CHAPTER  III 


A  bivariate  logit  is,  perhaps,  as  natural  to  consider  in  the  study 
of  a  mixture  of  two  stimulants  as  a  bivariate  normit,  even  though  very 
little  work  has  been  done  along  these  lines. 

Ashford  and  Smith  (39f,4)  ran  an  analysis  on  data  set  seven  assuming 
the  marginals  to  be  logistic.  They  fitted  the  data  to  a  model  using  a 
ccwplieated  response  region  without  explicitly  defining  the  mathematical 
model.  There  does  not  appear  to  have  been  any  other  examinations  of 
data  by  means  of  a  bivariate  logit  in  the  literature 

In  the  case  of  a  bivariate  logit,  the  first  consideration  is  the 
form  of  the  bivariate  distribution  to  be  used.  The  bivariate  logistic 
distribution  utilized  in  this  study  was 


![1  +  exp  {-  x)]  1[1  +  exp (-  y) ]  1 

•  {l  +  aQ[l  +  exp  (-  x)]  1[1  +  exp (-  y)] 

•  exp  (-  x-y)}  ,  — <  x  ,  y  <  00 

which  was  developed  by  Gumbel  (1961) .  The  density  function  is 


f3 (x,y) 


{exp (-  x  -  y) *[1  +  exp (-  x) ] 

_  a 

•  [1  +  exp (-  y) ]  }  •  {l  +  aQ[l  -  exp(-  x) 

-  exp (-  y)  +  exp (-  x  -  y)]/[l  +  exp(-  x) 

+  exp (-  y)  +  exp (-  x-y)]}  ,  -«  <  x  ,  y 
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The  correlation  coefficient  is 


(23) 


(24) 
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P  =  3a0/ir2  ,  (25) 

2 

where  -  1  <  a  <  1  :  thus  |p|  <  3/tt 

A  bivariate  logit  analysis  was  run  on  the  seven  data  sets  utilizing 
the  Gumbel  bivariate  logistic  distribution.  The  analytic  model  for  the 
bivariate  logit  analysis  was 

P(z1  ,  z  ,  0)  =  { 1  +  explB^Zj.  +  a1>]}  1  +  {l  +  exp[B2(z2  +  a2)]}  1 
"  {1  +  exp[B1(z1  +  a1)]}"1  {l  +  exp[B2(z2  +  a^]}"1 
■  +  aQ{l  +  exp(B1(z1  +  a1)]}  1  {l  +  exp[B2(z2  +  a2 ) ] }  1 

•  exptB^Zj^  +  a^)  +  B2(z2  +  a2)  ]  j  . 

,  z2  <  “  .  (26) 

A  resume  of  the  results  is  given  in  Table  2.  The  entries  of  Table  2  are 

the  same  as  those  of  Table  1. 


Data 

Set 

No. 

N 

No.  of 
Signi ficant 
Chi-squares 

r 

SSE  | 

.  1 

d.f . 

SSR 

d.f. 

R2 

1 

18 

5 

76.403  | 

13 

48773 

4 

.99844 

2 

10 

4 

29.262  | 

5 

574.90 

4 

.95157 

3 

24 

0 

15.603  | 

19 

3165.4 

4 

.99509 

4 

17 

2 

29.616  | 

12 

1365.2 

4 

.97877 

5 

12 

3 

30.169  | 

7 

675.69 

4 

.95724 

6 

15 

0 

11.976  | 

10 

20170 

4 

.99941 

7 

40 

2 

39.020  | 

35 

221.36 

4 

.85014 

TABLE  2 
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As  was  the  case  for  the  bivariate  normit,  tne  regression  was  found 
to  he  si  gni  f  i  r  anf  for-  each  of  the  seven  data  sets,  and  data  sets  three/ 
six/  and  seven  have  nonsignificant  chi-squares  indicating  no  significant 
departure  from  regression. 

For  each  data  set,  the  SSE  from  the  bivariate  logit  analysis  was 

larger  than  the  corresponding  SSE  from  the  bivariate  normit  analysis. 

2 

Similarly  R  from  the  bivariate  logit  analysis  for  each  data  set  was 

2 

smaller  than  the  corresponding  R  from  the  bivariate  normit  analysis. 

The  bivariate  logit  analysis  indicated  that  the  same  number  of  data  points 
differed  significantly  from  the  bivariate  logit  model  as  was  the  case  with 
bivariate  normit  model  for  each  data  set  with  the  exception  of  data  set 
two  (simple  similar  action) ,  and  data  set  five.  With  data  set  two,  the 
bivariate  logit  analysis  indicated  that  four  out  of  the  ten  data  points 
differed  significantly  from  the  bivariate  logit  model  as  compared  to  one 
out  of  ten  in  the  bivariate  normit  analysis.  With  data  set  five,  the 
bivariate  logit  analysis  indicated  that  three  out  of  the  twelve  data 
points  differed  significantly  from  the  bivariate  logit  model  as  compared 
to  two  out  of  twelve  in  the  bivariate  normit  analysis. 

On  the  whole,  the  bivariate  logit  analysis  did  not  do  as  well  as 
the  bivariate  normit  analysis,  although  it  did  nearly  as  well  with  six 
out  of  seven  of  the  data  sets.  It  would  seem  likely  that  the  main  reason 
that  the  bivariate  logit  model  did  not  do  as  well  was  due  to  the  fact 
that  the  correlation  coefficient  of  the  model  employed  was  restricted  so 
that  | p |  <  0.30396  ,  approximately,  and  |p|  >  0.30396  for  all  seven  data 
sets  in  the  bivariate  normit  analysis.  It  would  be  useful  to  extend  this 
investigation  to  include  a  bivariate  logit  model  where  the  correlation 
coefficient  is  not  so  restricted,  i.e.,  where  -  1  <  p  <  1  inclusive. 


CHAPTER  IV 


In  this  chapter  the  assumption  that  the  marginals  follow  the  Burr 
distribution  will  be  made.  This  is  a  somewhat  more  general  assumption 
than  the  assumption  that  the  marginals  are  normal  (or  logistic) . 

The  general  system  of  distributions  referred  to  here  was  first 
given  by  Burr  (1942).  Using  as  an  expression  for  the  distribution  function 

(  1  -  (1  +  xb)  p  x  >  0  ;  b  ,  p  >  0 
F(x)  =  <  (27) 

(  0  X  <  0 

F(x)  covers  an  important  region  of  the  standardized  third  and  fourth  central 

moments  in  the  following  sense.  Figure  1  shows  that  the  system  covers  a 

large  portion  of  the  curve-shape  characteristics  for  Types  I,  III,  IV,  and 

2 

VI  of  the  Pearson  system.  Figure  (1)  is  drawn  with  coordinates  »  0^ 

2  ^ 

and  6  =  (2a4  -  3a3  -  6)/(a4  +  3),  where  is  the  i  standardized  central 
moment.  The  regions  covered  by  the  Pearson  Types  I  (or  beta) ,  IV,  and  VI 
are  indicated,  as  well  as  Type  III  (or  gamma)  which  lies  on  a  curve,  and 
the  normal,  logistic, rectangular,  and  exponential  distributions  which  are 
represented  by  points.  The  subscript  B  refers  to  bell  shaped  functions 
and  J  to  J  shaped  functions.  It  can  be  seen  from  Figure  (1)  that  this 
system  of  distributions  is  quite  general. 

Takahasi  (1965)  developed  a  multivariate  Burr  distribution  by  using 
the  fact  that  a  Burr  distribution  is  a  compound  Weibull  distribution  with 
a  gamma-distribution  as  a  compounder.  That  is,  if 
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It  should  be  noted  at  this  point  that  the  are  equal  to  one  in  the  Burr 
distribution  as  given  by  Burr  (1942) .  If  x^  is  set  equal  to  B^(z^  +  <0  , 
it  is  easily  seen  that  the  r.'s  are  redundant.  In  addition,  if  the  b.'s 

i  l 

and  p  are  held  constant,  e.g. ,  the  third  and  fourth  standardized  central 
moments  can  be  set  equal  to  those  of  the  normal  distribution  by  proper 
choice  of  the  b^'s  and  p  ,  then  the  correlation  coefficient  is  a  constant. 

It  was  attempted  to  find  a  form  of  a  bivariate  Burr  distribution  such 
that  the  correlation  coefficient  would  not  be  a  fixed  constant.  The  form 


developed  by  the  author  is 
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0 


elsewhere 


(35) 
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The  correlation  coefficient  is 


v,  ■  r(p  -  ^)r(p  -  ^H1  ^>(i  ■ 

i1  -  £)  /i(- + 1  -  ^  - 1  *  r3  ”> + 1  - 1  -  *) 


2fi(p  •  *T  '  1  +  t  '•  P  +  1  }  1  -  * 


[r(1+^r(p-^)r,p>-r2(1+^r2(p-^)]- 
,r(l  +  ^)r(P  '  E|)r(P)  '  f2(‘  *  ^)rtE  "  £)]f  ' 


where  (a  ,B  ;y  ,z)  is  Gauss'  hypergeometric  function.  If  r  =  1  ,  then 


p  =  0  . 

X1X2 


A  bivariate  Burrit  analysis  was  run  on  the  seven  sets  of  data  using 

the  bivariate  Burr  distribution  described  above.  The  analytic  model  for 

the  bivariate  Burrit  analysis  was 

/  b  b 

P(z1  ,  z2  ,  0)  -  1  -  |  1  +  +  *1)1  1  +  [B2(z2  +  a2)]  2 

bl  b2)-D 

+  rtB.^^  +  a^)  ]  tB2(z?  +  a2^  |  •  "ai  i  2-1  <  “ 


“a2  -  z2  <  °° 


elsewhere  . 


A  resume  of  the  results  is  given  in  Table  3.  The  rows  corresponding  to 
eight  parameters  to  be  estimated  constitute  the  general  case  of  the  Burrit 
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analysis;  the  rows  corresponding  to  seven  parameters  to  be  estimated 

correspond  to  the  special  case  with  r  =  0  which  reduces  to  a  Burrit  analysis 

using  the  bivariate  Burr  distribution  developed  by  Takahasi  (1965).  The 

rows  corresponding  to  five  and  four  parameters  to  be  estimated  have  =  0  , 

=  3  (the  third  and  fourth  standardized  central  moments),  which  are  the 

same  as  the  normal  distributions'  a_  and  a.  .  The  first  of  these  is  a 

3  4 

special  case  of  the  general  Burrit  analysis  and  the  second,  a  special  case 
of  the  Burrit  analysis  using  the  Takahasi  bivariate  Burr  distribution. 

As  was  the  case  for  both  the  bivariate  normit  and  the  bivariate  logit, 
the  regression  was  found  to  be  significant  for  each  of  the  seven  data  sets 
for  all  of  the  bivariate  Burrit  analyses  (four  on  each  data  set) .  The 
chi-square  test  was  insignificant,  indicating  no  significant  departure 
from  regression  for  data  set  three  with  the  general  Burrit  analysis ,  for 
data  set  six  for  all  but  the  analysis  with  four  parameters  to  be  estimated, 
and  for  data  set  seven  for  all  four  of  the  analyses. 

The  SSE  from  the  general  case  of  the  bivariate  Burrit  analysis  was 
significantly  smaller  than  that  from  the  bivariate  normit  analysis  only 
with  the  synergistic  data  (data  set  one).  In  no  other  case  is  there  any 
indication  that  the  bivariate  Burrit  model  is  better  than  the  bivariate 
normit  model  in  the  actual  fitting  of  these  data  to  a  model. 

Each  SSE  from  the  bivariate  Burrit  analyses  utilizing  the  bivariate 
Burr  developed  in  this  paper  is  significantly  smaller  than  the  corresponding 
SSE  from  the  analyses  utilizing  the  Takahasi  bivariate  Burr  distribution 
in  all  but  three  cases:  both  cases  with  data  set  seven,  and  the  first  case 
with  data  set  five  (the  case  corresponding  to  the  two  analyses  with  eight 
and  seven  parameters  to  be  estimated) .  On  the  basis  of  these  analyses  it 
would  seem  that  the  bivariate  Burr  developed  in  this  paper  would  be,  in 
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general,  more  useful  in  application  than  the  form  developed  by  Takahasi. 

T-—  -  Jj;  A.  1  £  a. ~  _  x-  -  J  i.1.  .i.  Xt  -  — - J 1  - - JL..XJ  £-  — 

XU  UUUX  <.Xt/M  X  Ir  iiiujf  •*-»'—  X«t^*  <~4  WilWW  WUW  tltUX^XiiUX  UlaltXXVWbXVUg  XWX 

synergistic  data,  as  characterized  by  and  ,  do  not  lie  in  the  same 
Pearson  curve  area  (see  figure  2) .  The  marginals  for  data  set  three  also 
display  this  characteristic  but  not  to  as  high  a  degree.  The  marginals 
for  data  sets  four  through  seven  are  all  clustered  around  the  normal  dis¬ 
tribution.  The  fact  that  the  assumption  that  the  marginals  are  Burr  dis¬ 
tribution  does  allow  given  marginal  to  have  curve  shape  characteristics 
different  from  that  of  the  other  marginal  suggests  that  the  bivariate 
Burrit  analysis  may  be  well  adapted  for  the  analysis  of  data  where  the 
marginal  distributions  do  not  belong  to  the  same  family,  e.g.  the  family 
of  normal  distributions. 

In  summary,  the  bivariate  analyses  utilizing  the  general  form  indicated 
by  equation  (16)  seem  to  do  quite  well  with  a  diversity  of  mixtures  of 
stimulants  as  is  evidenced  by  the  seven  sets  of  data  which  have  been  ana¬ 
lyzed  in  this  paper.  The  bivariate  nonnit  model  and  the  bivariate  Burrit 
model  (general  case,  i.e.,  the  case  with  eig.»t  parameters  to  be  estimated) 
seem  to  be  best  suited  for  these  types  of  analyses.  The  bivariate  normit 
model  would  have  to  be  recommended  if  the  number  of  parameters  to  be  esti¬ 
mated  is  of  concern,  but  otherwise  the  bivariate  Burrit  model  could  well 
be  the  best  model  for  these  types  of  analyses. 
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Figure  2.  Expanded  portion  of  the  coverage  8^  ,  6  space 
The  x's  mark  six  of  the  sample  population  points  (8^  , 
from  the  data  sets  analyzed  in  this  paper.  N^(N=1,  2,  3 
i=a,  b)  refers  to  the  i^1  marginal  of  the  data  set. 


APPENDIX  I 


Data  of  Kagy  and  Richardson  (1936) :  The  combined  action  of  2-4 
dinitro-6-cyclohexylphenol  and  petroleum  oil  sprayed  in  emulsions  against 
eggs  of  a  plant  bug  (Lygaeus  Kalmii  Stal) .  Th<  data  as  described  by  Kagy 
and  Richardson,  the  translated  data,  and  the  analyses  on  this  set  of  data 
(data  set  one)  are  in  this  appendix. 
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Here  (N^ (it*1  Net  kill  %)/100)  was  rounded  off  to  the  nearest  integer  — 
wh:.i..i  should  be  the  number  that  responded  to  the  i^  mixture  of  stim¬ 
ulants,  and  then  p^  was  computed  as  rVN^. 


I 
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BIVARIATE  NORMIT  ANALYSIS 


Parameter  Estimates 
a  =  8.640 

B  =  0.935 

a2  =  5.583 

B2  =  1.489 

p  =  -0.379 


Chi-square  Analysis  Table 


1 . — - 

— 

1 

source 

d.  f . 

Due  to  Model 

SSR 

4 

66041 

Departure  from  Model 

SSE 

13 

67.029 

TOTAL 

SST 

17 

66107 

P. 

1 

.0667 

.1016 

.4008 

.4049 

.5866 

.6417 

.1000 

.1180 

.5971 

.4637 

.7223 

.7231 

.3021 

.3596 

.7370 

.7675 

.9040 

.9277 

.5868 

.5859 

.9401 

.8986 

.9722 

.9768 

.8125 

.7986 

.9714 

.9691 

.9965 

.9950 

.8681 

.8865 

.9948 

.9870 

.9667 

.9536 

Coefficient  of  Determination  R  =  .99899 

Residual  Analysis 
P.  Residual  Chi-square 


-.0349 

-.0041 

-.0551 

-.0180 

.1334 

-.0008 

-.0575 

-.0305 

-.0237 

-.0009 

.0415 

-.0046 

-.0139 

.0023 

.0015 

-.0184 

.0078 

.0131 


3.201 
0.033 
6.317 
0.750 
34.266 
0.002 
4.139 
2.490 
4.024 
0.001 
7.266 
0.274 
0.  346 
0.069 
0.137 
0.973 
1.809 
0.931 


<  ed 
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BIVARIATE  LOGIT  ANALYSIS 


Parameter  Estimates 
aQ  =  -0.938 
=  9.220 

B  =  -1.690 
a2  =  3.762 

B2  =  -2.413 


Chi-square  Analysis  Table 


source 

d.f . 

Due  to  Model 

SSR 

4 

48773 

Departure  from  Model 

SSE 

13 

76.402 

TOTAL 

SST 

17 

48849 

2 

Coefficient  of  Determination  R  =  .99844 


Residual  Analysis 


Pi 

P. 

1 

Residual 

Chi-square 

.0667 

.1155 

-.0488 

5.597 

.4008 

.4103 

-.0095 

0.177 

.5866 

.6492 

-.0626 

8.246 

.1000 

.1352 

-.0352 

2.540 

.5971 

.4601 

.1370 

36.200 

.7223 

.7120 

.0103 

0.247 

.3021 

.3438 

-.0417 

2.220 

.7370 

.7627 

-.0257 

1.749 

.9040 

.9303 

-.0263 

5.117 

.5868 

.5837 

.0031 

0.011 

.9401 

.9013 

.0388 

6.503 

.9722 

.9774 

-.0052 

0.348 

.8125 

.8087 

.0038 

0.027 

.9714 

.9658 

.0056 

0.365 

.9965 

.9929 

.0036 

0.520 

.8681 

.8917 

-.0236 

1.654 

.9948 

.9821 

.0127 

3.527 

.9667 

.9504 

.0163 

1.354 
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BIVARIATE  BURRIT  ANALYSES 

1:  General  Case  -  Eight  Parameters  to  De  Estimated 


Parameter  Estimates 
r  =  3.556 
b1  =  9.643 
b2  =  1.773 
p  =  4.813 
aL  =  18.073 
B  =  0.094 
a2  =  4.877 
B2  =  0.318 


Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

7 

67964 

Departure  from  Model 

SSE 

10 

41.646 

TOTAL 

SST 

17 

68006 

Coefficient  of  Determination 


.99939 


Residual  Analysis 


Pi 

P. 

1 

Residual 

Chi-square 

.0667 

.0602 

.0065 

0.181 

.4008 

.4273 

-.0265 

1.377 

.5866 

.6289 

-.0423 

3.676 

.1000 

.1018 

-.0018 

0.008 

.5971 

.5054 

.0917 

16.107 

.7223 

.7202 

.0021 

0.010 

.3021 

.3471 

-.0450 

2.575 

.7370 

.7743 

-.0373 

3.821 

.9040 

.9211 

-.0171 

1.935 

.5868 

.5627 

.0241 

0.682 

.9401 

.9010 

.0390 

6.580 

.9722 

.9755 

-.0033 

0.130 

.8125 

.7823 

.0302 

1.540 

.9714 

.9711 

.0003 

0.014 

.9965 

.9950 

.0015 

0.137 

.8681 

.8781 

-.0100 

0.271 

.9948 

.9883 

.0066 

1.420 

.9667 

.9516 

.0151 

1.194 
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2:  Takahasi  Burr  -  r  =  0  ;  Seven  Parameters  to  be  Estimated 


Parameter  Estimates 


16.156 

0.113 

4.869 

0.291 


Chi-square  Analysis  Table 


source 


Due  to  model 
Departure  from  Model 
TOTAL 


Coefficient  of  Determination 


45 
64.694 
45 


Residual  Analysis 


P. 

r 

Residual 

Chi-square 

.0667 

.0592 

.0075 

0.245 

.4008 

.4480 

-.0472 

4.318 

.5866 

.6599 

-.0733 

11.463 

.1000 

.0921 

.0079 

0.181 

.5971 

.5001 

.0970 

18.015 

.7223 

.7115 

.0108 

0.272 

.3021 

.3552 

-.0531 

3.544 

.7370 

.7466 

-.0096 

0.231 

.9040 

.8874 

.0166 

1.319 

.5868 

.5958 

-.0090 

0.097 

.9401 

.8844 

.0557 

0.001 

.9722 

.9581 

.0141 

1.420 

.8125 

.8261 

-.0136 

0.370 

.9714 

.9670 

.0044 

0.230 

.9965 

.9904 

.0061 

1.111 

.8681 

.9151 

-.0469 

8.167 

.9948 

.9875 

.0073 

1.656 
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3 1  ct3  =  0  ;  a4  =  3  (Third  and  Fourth  Standardized  Central  Moments)  t 
b  =  b2  =  4.874  ;  p  =  6.158  ;  Five  Parameters  to  be  Estimated 


Parameter  Estimates 
r  =  4.383 

a  =  13.485 
=  0.153 

a2  =  6.426 

B2  =  0.242 


Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

4 

70566 

Departure  from  Model 

SSE 

13 

70.017 

TOTAL 

SST 

17 

70636 

2 

Coefficient  of  Determination  R  =  .99901 


Residual  Analysis 


Pi 

P. 

1 

Residual 

Chi-square 

.0667 

.1057 

-.0390 

3.861 

.4008 

.4058 

-.0050 

0.049 

.  5866 

.6443 

-.0577 

6.967 

.1000 

.1221 

-.0221 

1.090 

.5971 

.4603 

.1358 

36.090 

.7223 

.7168 

.0056 

0.073 

.3021 

.3594 

-.0573 

4.107 

,7370 

.7592 

-.0222 

1.296 

.9040 

.9271 

-.0231 

3.794 

.  5368 

.5829 

.0039 

1.873 

.9401 

.8964 

.0437 

7.906 

.9722 

.9782 

-.0060 

0.483 

.8125 

.7985 

.0140 

0.353 

.9714 

.9691 

.0023 

0.068 

.9965 

.9954 

.0011 

0.073 

.8681 

.8876 

-.0195 

1.100 

.9948 

.9869 

.0079 

1.847 

.9667 

.9543 

.0124 

0.843 
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4:  Takahasi  Burr  -  r  =  0  ;  =  0  ;  =  3  (Third  and  Fourth 

Standardized  Central  Moments);  b^  =  =  4.874  ; 

p  =  6.158  ;  Four  Parameters  to  be  Estimated 

Parameter  Estimates 
ax  =  13.460 
B  =  0.160 

a2  =  6.265 

B2  =  0.259 

2 

Coefficient  of  Determination  R  =  .99733 


Residual  Analysis 


pi 

P. 

1 

Residual 

Chi-square 

.0667 

.0957 

-.0290 

2.333 

.4008 

.4125 

-.0117 

0.270 

.5866 

.6674 

-.0808 

14.102 

.1000 

.1140 

-.0140 

0.463 

.5971 

.4596 

.1375 

36,464 

.7223 

.7144 

.0079 

0.147 

.3021 

.3799 

-.0778 

7.395 

.7370 

.7378 

-.0008 

14.954 

.9040 

.8911 

.0129 

0.822 

.5868 

.6188 

-.0319 

1.246 

.9401 

.8789 

.0612 

13.499 

.9722 

.9564 

.0158 

1.723 

.8125 

.8321 

-.0196 

0.796 

.9714 

.9611 

.0103 

1.086 

.9965 

.9877 

.0088 

1.844 

.8681 

.9128 

-.0446 

7.208 

.9948 

.9830 

.0118 

3.208 

.9667 

.9680 

.0013 

0.013 

Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

3 

34582 

Departure  from  Model 

SSE 

14 

92.619 

TOTAL 

SST 

17 

34674 

APPENDIX  II 


Data  of  Plackett  and  Hewlett  (1952)  :  The  toxity  to  Tribolium 
castaneum  of  D.D.T. ,  methoxychlor  (MOC) ,  and  combinations  of  the  two 
applied  in  Shell  Oil  P31  as  films  on  filter  paper,  six-day  exposures. 
The  data  as  described  by  Plackett  and  Hewlett,  the  translated  data,  and 
the  analyses  on  this  set  of  data  (data  set  two)  are  in  this  appendix. 
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DATA  AS  DESCRIBED  BY  PLACKETT  AND  HEWLET' 


TRANSLATED  DATA 


BIVARIATE  NORM IT  ANALYST I 


Parameter  Estimates 


Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

4 

598.86 

Departure  from  Model 

SSE 

5 

21.775 

TOTAL 

SST 

9 

620.63 

0.9999 


Coefficient  of  Determination 


Residual  Analysis 


Pi 

P. 

1 

'  Residual 

Chi-square 

.0754 

.0844 

-.0090 

0.210 

.2973 

.3693 

-.0720 

3.292 

.7789 

.7606 

.0183 

0.365 

.1450 

.1920 

.0470 

2.848 

.2600 

.2764 

-.0164 

0.203 

,6358 

.5613 

.0745 

3.405 

.4362 

.4491 

-.0129 

0.100 

.6622 

.  5335 

.1287 

9.845 

.7867 

.8184 

-.0317 

1.013 

.7085 

.7306 

-.0221 

0.494 

Parameter  Estimates 
aQ  =  -1.000 
a  5.448 
B1  -  -1.776 
a2  =  4.645 

B2  =  -2.559 


.0754 

.2973 

.7789 

.1450 

.2600 

.6358 

.4362 

.6622 

.7867 

.7085 
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5IVARIATE  LOGIT  ANALYSIS 


Chi-Square  Analysis  Table 


source 

. 

d.f. 

Due  to  Model 

SSR 

4 

574.90 

Departure  from  Model 

SSE 

5 

29.262 

TOTAL 

SST 

9 

604.16 

2 

Coefficient  of  Determination  R  =  .95157 


Residual  Analysis 


P. 

1 

Residual 

Chi-square 

.0958 

-.0204 

0.960 

.3846 

-.0873 

4.760 

.7865 

-.0076 

0.068 

.2041 

-.0591 

4.297 

.2944 

-.0344 

0.856 

.5486 

.0872 

4.638 

.4675 

-.0313 

0.586 

.5401 

.1221 

8.811 

.7312 

.0555 

2.351 

.7504 

-.0419 

1.864 

w> 


BIVARIATE  BURRIT  ANALYSES 

1:  General  Case  -  Eight  Parameters  to  be  Estimated 

Parameter  Estimates 
r  =  5.271 
bx  =  2.351 
b2  =  6.755 
p  =  4.271 
a  =  7.252 
B  =  0.270 
a^  =  8.206 

2  "  °-216 

Residual  Analysis 


1 

i 

i 

i 

i 
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2:  Takahasi  Burr  -  r  =  0  ;  Seven  Parameters  to  be  Estimated 


Parameter  Estimates 


b  =  0.961 


b2  =  5.064 


p  =  3.265 


ax  =  6. 399 


B  =  0.298 


a2  =  7.263 


B2  =  0.290 


source 

d.f. 

Due  to  Model 

SSR 

6 

544.20 

Departure  from  Model 

SSE 

3 

38.249 

TOTAL 

SST 

9 

582.45 

Coefficient  of  Determination 


R  =  .93433 


Residual  Analysis 


fi 

P. 

l 

Residual 

Chi-square 

0754 

.0987 

-.0233 

1.217 

2973 

.4096 

-.1123 

7.723 

7789 

.7928 

-.0139 

0.234 

1450 

.1774 

-.0324 

1.440 

2600 

.2542 

.0058 

0.026 

6358 

.5004 

.1354 

11.071 

4362 

.5490 

-.1128 

7.662 

6622 

.5844 

.0778 

3.685 

7867 

.7037 

.0830 

4.953 

7085 

.7240 

-.0155 

0.239 

<11  > 
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3:  ct  =0;a  =3  (Third  and  Fourth  standardized  Central  Moments)  > 

3  4 

b1  =  fc>2  =  4.874  ;  p  =  6.158  j  Five  Parameters  to  be  Estimated 


Parameter  Estimates 
r  =  7.158 
=  9.108 
B  =  0.176 
a2  =  7.272 
B2  =  0.245 


Chi-square  Analysis  Table 


source 

d.f . 

Due  to  Model 

SSR 

4 

577.20 

Departure  from  Model 

SSE 

5 

26.296 

TOTAL 

SST 

9 

603.50 

2 

Coefficient  of  Determination  R  -  .95643 


Residual  Analysis 


Pi 

Pi 

Residual 

Chi-square 

.0754 

.0946 

-.0192 

0.854 

.2973 

.3869 

-.0896 

5.005 

.7789 

.7805 

-.0016 

0.003 

.1450 

.2018 

-.0568 

4.002 

.2600 

.2928 

-.0328 

0,780 

.6358 

.5586 

.0772 

3.653 

.4362 

.4628 

-.0266 

0.424 

.6622 

.5410 

.1212 

8.757 

.7867 

.7490 

.0377 

1.133 

.7085 

.7484 

-.0399 

1.685 

w> 
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4:  Takahasi  Burr  -  r  =  0  ;  «3  =  0  ;  ct^  -  3  (Third  and  Fourth 
Standardized  Central  Moments);  =  k>2  =  4.874  ; 
p  =  6.158  ;  Four  Parameters  to  be  Estimated 


Parameter  Estimates 


a 

B 

a 


1 

1 

2 

2 


9.196 

0.174 

7.270 

0.248 


Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

3 

559.09 

Departure  from  Model 

SSE 

6 

46.537 

TOTAL 

SST 

9 

605.63 

2 

Coefficient  of  Determination  R  =  .92316 


Residual  Analysis 


Pi 

P. 

1 

Residual 

Chi-square 

.0754 

.0990 

-.0236 

1.239 

.2973 

.4013 

-.1040 

6.666 

.7789 

.7952 

-.0163 

0.324 

.1450 

.2165 

-.0715 

6.027 

.  2600 

.2912 

-.0312 

0.707 

.6358 

.5219 

.1139 

7.655 

.4362 

.4798 

-.0436 

1.135 

.6622 

.5264 

.1358 

10.953 

.7867 

.6727 

.1140 

8.857 

.7085 

.7590 

-.0505 

_ 

2.775 

APPENDIX  III 


Data  of  Hewlett  and  Plackett  (1950) :  A  study  of  six  day  toxicity  to 
beetles  (Tribolium  castaneum)  of  direct  sprays  of  Pyrethins,  D.D.T.,  and 
the  two  together  in  Shell  Oil  P31.  The  data  as  reproduced  by  Zeigler  and 
Moore  (1966)  ,  and  the  analyses  on  this  set  of  data  (data  set  three)  are  in 
this  appendix. 
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DATA  AS  REPRODUCED  BY  ZEIGLER  AND  MOORE 


Insecticide  (mg./10  sq.  cm.) 


1.2%  w/v 

2.52 

Pyrethins 

3.30 

4.25 

5.33 

7.15 

9.53 

12.28 

15.58 

2.0%  w/v 

2.45 

D.D.T. 

3.18 

4.25 

5.48 

7.24 

9.54 

12.36 

15.54 

1.2%  w/v 
Pyre thins 


2.0%  w/v 
D.D.T. 


2.74 
3.20 
4.  10 
5.34 
7.11 
9.60 
12.45 


.03024 

.03960 

.05100 

.06396 

.08580 

.11436 

.14739 

.18696 


.02964 

.03840 

.04920 

.06408 

.08532 

.11520 

.14940 


49 

.1633 

50 

.1600 

50 

.3200 

50 

.4200 

50 

.5000 

50 

.5600 

50 

.7000 

50 

.7400 

50 

.2800 

49 

.3673 

50 

.4400 

50 

.7200 

50 

.8400 

50 

.9000 

15.65 


18780 


3130 


5i 


<  Q. 
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BIVARIATE  NORMIT  ANALYSIS 


Parameter  Estimates 
a^  =  2.827 

Bx  =  1.231 

a2  =  1.698 

B2  =  0.882 

=  -0.686 


Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

4 

11176 

Departure  from 

Model 

SSE 

19 

11.805 

TOTAL 

SST 

23 

11188 

fi 

P. 

1 

.0625 

.0696 

.0625 

.1257 

.1800 

.2017 

.3200 

.2888 

.4000 

.4225 

.6000 

.5629 

.7551 

.6809 

.7000 

.7773 

.1633 

.1677 

.1600 

.2318 

.3200 

.3166 

.4200 

.4002 

.5000 

.4972 

.5600 

.5934 

.7000 

.6790 

.  7400 

.7476 

.2800 

.2358 

.3673 

.3506 

.4400 

.4896 

.7200 

.6559 

.8400 

.8215 

.9000 

.9364 

1.0000 

.9816 

1.0000 

.9954 

Coefficient,  of  Determination  R  =  .99894 

Resi dual  Analysis 

Residual  Chi-square 
696  -.0071  0.037 


-.0632 

-.0217 

.0312 

-.0225 

.0371 

.0742 

-.0773 

-.0044 

-.0718 

.0034 

.0198 

.0028 

-.0334 

.0210 

-.0076 

.0442 

.0167 

-.0496 

.0642 

.0185 

-.0364 

.0184 

.0046 


0.037 

1.744 

0.146 

0.237 

0.104 

0.280 

1.241 

1.727 

0.007 

1.447 

0.003 

0.082 

0.002 

0.231 

0.101 

0.015 

0.542 

0.060 

0.492 

0.912 

0.117 

1.109 

0.939 

0.231 
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BIVARIATE  LOGIT  ANALYSIS 


Parameter  Estimates 
aQ  -  -1.000 
a  =  2.338 

B  =  -2.090 
a2  =  1 . 960 

§2  =  -1.520 


Chi-squarc  Analysis  Tabic 


source 

d.f. 

Due  to  Model 

SSR 

4 

3165.4 

Departure  from  Model 

SSE 

19 

15.603 

TOTAL 

SST 

23 

3181.0 

2 

Coefficient  of  Determination  R  =  .99509 


Residual  Analysis 


pi 

P. 

1 

Residual 

Chi-square 

.0625 

.0813 

-.0188 

0.227 

.0625 

.1346 

-.0720 

2.140 

.  1800 

.2087 

-.0287 

0.250 

.  3200 

.2975 

.0225 

0.121 

.4000 

.4390 

-.0390 

0.309 

.6000 

.5879 

.0121 

0.030 

.7551 

.7079 

.0472 

0.528 

.7000 

.7994 

-.0994 

3.082. 

.1633 

.1674 

-.0041 

0.006 

.  1600 

.2301 

-.0701 

1.385 

.  3200 

.  3171 

.0029 

0.002 

.4200 

.4059 

.0141 

0.041 

.5000 

.5105 

-.0105 

0.022 

.5600 

.6133 

-.0533 

0.600 

.  7000 

.7016 

-.0016 

0.001 

.7400 

.7690 

-.0290 

0.237 

.2800 

.2442 

.0358 

0.347 

.  3673 

.  3493 

.0180 

0.070 

.4400 

.4754 

-.0354 

0.252 

.7200 

.6266 

.0934 

1.865 

.8400 

.  7816 

.0584 

0.999 

.9000 

.9002 

-.0002 

0.  000 

1.0000 

.9578 

.0422 

2.201 

1.0000 

.9826 

.0174 

0.888 

— 

JU> 
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BIVARIATE  BURRIT  ANALYSES 


uelieiai 


K/%  Pof  rna+- 10 rl 


Parameter  Estimates 


r  =  6.323 
6  =  3.933 

b2  =  2.941 
p  =  5.323 
a  =  4.863 
B  =  0.239 
2  =  4.760 
B2  =  0.183 


Chi-sauare  Analysis  Table 


source 

d.f . 

Due  to  Model 

SSR 

7 

4258.4 

Departure  from  Model 

SSE 

16 

13.227 

TOTAL 

SST 

23 

4271.6 

Coefficient  of  Determination 


R2  =  .99690 


Residual  Analysis 


Pi 

P, 

1 

Residual 

Chi-square 

.0625 

.0622 

.0003 

0.000 

.0625 

.1216 

-.0591 

1.569 

.1800 

.2025 

-.0225 

0.156 

.3200 

.2947 

.0253 

0.154 

.4000 

.4346 

-.0346 

0.244 

.6000 

.5784 

.0216 

0.095 

.7551 

.6963 

.0588 

0.802 

.7000 

.7898 

-.0898 

2.429 

.1633 

.  1666 

-.0033 

0.004 

.1600 

.2384 

-.0784 

1.691 

.  3200 

.  3298 

-.0098 

0.022 

.4200 

.4161 

,0039 

0.003 

.5000 

.5120 

-.0120 

0.029 

.5600 

.6034 

-.0434 

0.393 

.7000 

.6820 

.0180 

0.075 

.7400 

.7438 

-.0038 

0.0 

.2800 

.2261 

.0540 

0.832 

.3673 

.3475 

.0198 

0.085 

.4400 

.4875 

-.0475 

0.451 

.7200 

.6463 

.0737 

J  .187 

.8400 

.8007 

.0393 

0.484 

.9000 

.9135 

-.0135 

0.115 

1.0000 

.9660 

.0340 

1.759 

1.0000 

.9873 

.0127 

0.644 

i 


2:  Takahasi  Burr 


0  ;  Seven  Parameters  to  be  Estimated 


Parameter  Estimates 


a.  a  ji 


t>2  =  2.939 


p  =  6.740 


a^  =  4.675 


Bx  =  0.243 


a2  =  4.701 


B2  =  0.176 


Chi-square  Anai ys is ^ Table 


source 

d.f . 

Due  to  Model 

SSR 

& 

1178.8 

Departure  trom  Model 

SSE 

17 

29.737 

TOTAL 

SST 

23 

1208.5 

Coefficient  of  Determination  R  =  .97539 


Residual  Analysis 


?i 

P, 

1 

Residual 

Chi-square 

.0625 

.0540 

.0085 

0.067 

.0625 

.1148 

-.0523 

1.292 

.1800 

.2011 

-.0211 

0.138 

.  3200 

.3017 

.0184 

0.080 

.4000 

.4552 

-.0552 

0.614 

.6000 

.6110 

-.0109 

0.025 

.7551 

.7346 

.0205 

0.106 

.7000 

.8282 

-.1282 

5.772 

.1633 

.1707 

-.0074 

0.019 

.1600 

.2469 

-.0869 

2.032 

.3200 

.3444 

-.0244 

0.132 

.4200 

.4361 

-.0161 

0.053 

.5000 

.5373 

-.0373 

0.280 

.5600 

.6326 

-.0725 

1.132 

.7000 

.7132 

-.0132 

0.042 

.7400 

.7752 

-.0352 

0.355 

.2800 

.2137 

.  (;  j63 

1.310 

.  3673 

.3256 

.0417 

0.387 

.4400 

.4503 

-.0103 

0.022 

.  7200 

.5888 

.1312 

3.554 

.8400 

.7267 

.1134 

3.234 

.9000 

.8406 

.0594 

1.318 

1.0000 

.9088 

.0912 

5.019 

1.0000 

.9478 

.0522 

2.754 

m 

I 


I 

! 


i 


I 

| 

i 
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3:  -i  =  0  ;  i  3  (Third  and  Fourth  Standardized  Central  Moments)  ; 

3  4 

b,  =  b  =  4.874  ;  i-  -  0.158  ;  Five  Parameters  to  be  Estimated 
1  2 

Parameter  Estimates 
r  =  7.129 
ax  =  5.499 
Bx  =  0.200 
a  =  5.056 
§2  =  0.207 


Residual  Analysis 


pi 

P  , 

1 

Residual 

Chi-square 

.0625 

.0679 

-.0054 

0.022 

.0625 

.1215 

-.0590 

1.566 

.1800 

.1937 

-.0137 

0.060 

.3200 

.2771 

.0429 

0.460 

.4000 

.4076 

-.0076 

0.012 

.6000 

.5483 

.0517 

0.540 

.7551 

.6695 

.0856 

1.624 

.7000 

.7699 

-.0699 

1.378 

.1633 

.0884 

.0749 

3.408 

.1600 

.1525 

.0075 

0.022 

.3200 

.2532 

.0668 

1.178 

.4200 

.3654 

.0546 

0.644 

.5000 

.5046 

-.0046 

0.004 

.5600 

.6438 

-.0338 

1.530 

.7000 

.7602 

-.0602 

0.995 

.7400 

.8431 

-.1031 

4.014 

.2800 

.1542 

.1258 

6.062 

.  3673 

.2656 

.1018 

2.601 

.4400 

.4094 

.0306 

0.194 

.7200 

.5903 

.1297 

3.477 

.8400 

.7803 

.0597 

1.041 

.9000 

.9197 

-.0197 

0.263 

1.0000 

.9767 

.0233 

1.193 

1.0000 

.9942 

.0058 

0.291 

Chi-square  Analysis  Table 


source 

— 

d.f. 

Due  to  Model 

S9R 

4 

8668.4 

Departure  from  Model 

SSE 

19 

32.579 

TOTAL 

SST 

23 

8701.0 

2 

Coefficient  of  Determination  R  =  .99626 


<1:  Takahar.i  Hurr  -  r  -  0  ;  ,i  ,  -  0  ;  i,  =  3  (Third  and  Fourth 

J  4 

Standardized  Central  Moments);  ~  4.874  ; 

p  =  6.158  ;  Four  Parameters  to  be  Estimated 

Parameter  Estimates 
ax  =  4.932 
B  =  0.246 
a2  =  4.585 
B2  =  0.237 


Residual  Analysis 


pi 

P. 

1 

Residual 

Chi-square 

.0625 

.0377 

.0248 

0.813 

.0625 

.0849 

-.0224 

0.309 

.  1800 

.1588 

.0212 

0.168 

.  3200 

.2533 

.0667 

1.176 

.4000 

.4114 

-.0114 

0.027 

.6000 

.5849 

.0151 

0.047 

.  7551 

.7280 

.0271 

0.182 

.7000 

.8357 

-.1357 

6.706 

.1633 

.0482 

.1151 

14.133 

.1600 

.0989 

.0611 

2.094 

.3200 

.1906 

.1294 

5.423 

.4200 

.304  3 

.1157 

3.162 

.5000 

.4569 

.0431 

0.374 

.5600 

.6174 

-.0573 

0.696 

.7000 

.7533 

-.0533 

0.764 

,7400 

.8480 

-.1080 

4.521 

.2800 

.0827 

.1973 

25.672 

.3673 

.1693 

.1978 

13.626 

.4400 

.2940 

.1460 

5.132 

.7200 

.4632 

.2568 

13.259 

.8400 

.6569 

.1831 

7.440 

.9000 

.8241 

.0759 

1.987 

1.0000 

.9172 

.0828 

4.513 

1.0000 

.9625 

.0375 

1.947 

Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

3 

1309.0 

Departure  from  Model 

SSE 

20 

114.16 

TOTAL 

SST 

23 

1423.2 

2 

Coefficient  of  Detemination  R  =  .91978 


APPENDIX  IV 


DaLcj  uT  J.  T.  Mai: tin  (1042):  The  toxiciticc  to  Macros iphoniella 
sanborni  of  rotennne,  a  deguelin  concentrate,  and  of  a  mixture-  Tests 
of  17  November  1938.  Fivefold  replication.  Results  one  day  after  spraying. 
Medium  0.5%  saponin,  containing  alcohol.  Tattersfield  apparatus.  The 
data  as  described  by  Martin,  the  translated  data,  and  the  analyses  of  this 
set  of  data  (data  set  four)  are  in  this  appendix. 
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5b 


DATA  A3  DESCRIBED  BY  MARTIN 


CONCENTRATIONS  (mg./i.) 

Ni 

Number  of 

Insects 

Used 

Percent 

Mortality 

X(l) 

Rotenone 

X  (2) 

Deguelin 

Concentrate 

10.2 

0.0 

50 

88.0 

0.0 

49 

85.7 

I'. 

0.0 

46 

52.2 

3.8 

0.0 

48 

33.3 

2.6 

0.0 

50 

12.0 

0.0 

50.5 

48 

100.0 

0.0 

40.4 

50 

94.0 

0.0 

30.3 

49 

95.9 

0.0 

20.2 

48 

70.8 

0.0 

10.1 

48 

37.5 

0.0 

5.1 

49 

32.6 

5.1 

20.3 

50 

96.0 

4.0 

16.3 

46 

93.5 

3.0 

12.2 

48 

79.2 

2.0 

8.1 

46 

58.7 

1.0 

4.1 

46 

47.8 

0.5 

2.0 

47 

14.9 

TRANSLATED  DATA 


Pi 

.8800 

.8571 

.5217 

.3333 

.1200 

1.0000 

.9400 

.9592 

.7083 

.3750 

.3265 

.9600 

.9348 

.7917 

.5370 

.4783 

.1489 
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BIVARIATE  NORMIT  ANALYSIS 


Parameter  Estimates 
a  =  -2.775 
B  =  1.762 

a  -  -1.645 
B2  =  0.823 

(5  =  -0.530 


Chi-square  Analysis  Table 


source 

Due  to  Model 

1 

SSR 

d.f. 

4 

1656.0 

Departure  from  Model 

SSE 

12 

27.146 

TOTAL 

SST 

16 

1683.1 

2 

Coefficient  of  Determination  R  =  .98387 


Residual  Analysis 


Pi 

Pi 

Residual 

Chi-square 

.8800 

.9059 

-.0259 

0.392 

.8571 

.7940 

.0631 

1.193 

.5217 

.5377 

-.0160 

0.047 

.3333 

.  3359 

-.0026 

0.001 

.1200 

.  1374 

-.0174 

0.128 

1.0000 

.9432 

.0568 

2.893 

.9400 

.9190 

.0210 

0.297 

.9592 

.6773 

.0819 

3.053 

.7083 

.7962 

-.0879 

2.284 

.  3750 

.6017 

-.2267 

10.297 

.  3265 

.  3805 

-.0540 

0.605 

.9600 

.9649 

-.0049 

0.036 

.9348 

.9084 

.0264 

0.386 

.7917 

,7893 

.0024 

0.002 

.5870 

.5826 

.0044 

0.004 

.4783 

.3170 

.1613 

5.528 

.1489 

.  1506 

-.0017 

0.001 
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BIVARIATE  LOGIT  ANALYSIS 


Chi -souare  Analvsis 

Table 

a.  =  -1.000 

source 

d.f. 

a,  =  -1.568 

Due  to  Model  3SR 

4 

1365.2 

B  -  -  2  9£4 

Departure  from  Model  SSE 

12 

29.616 

TOTAL  S3T 

16 

1394.9 

a.  =  -2.029 

B2  =  -1.450 


Coefficient  of  Determination  R  =  .97877 


Residual  Analysis 


£i 

P, 

1 

Residual 

Chi-square 

.8800 

.9034 

-.0234 

0.315 

.8571 

.8026 

.0545 

0.918 

.5217 

.5453 

-.0236 

0.103 

.3333 

.3392 

-.0006 

0.000 

.1200 

.1400 

-.0200 

0.166 

1.0000 

.9397 

.0603 

3.082 

.9400 

.9185 

.0215 

0.309 

.9592 

.8813 

.0779 

2.84'’ 

.7083 

.8048 

-.0965 

2.848 

.3750 

.6014 

-.2264 

10.267 

.3265 

.3590 

-.0325 

0.225 

.9600 

.9505 

.0095 

0.095 

.9348 

.8865 

.0484 

1.069 

.7917 

.7672 

.0245 

0.161 

.5870 

.5722 

.0148 

0.041 

.4783 

.2985 

.1798 

7.099 

.1489 

.1352 

.0137 

0.076 
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BIVARIATE  BURRIT  ANALYSES 

1:  General  Case  -  Eight  Parameters  to  be  Estimated 

Parameter  Estimates 
r  =  7.109 
b  =  4.891 
b  =  4.511 
p  =6.109 
a^  =  0.695 
B  =  0.286 
a2  =  2.730 
B2  =  0.136 


Residual  Analysis 


pi 

P. 

1 

Residual 

Chi-square 

.8800 

.9125 

-.0325 

0.660 

.8571 

.8021 

.0550 

0.934 

.5217 

.5430 

-.0213 

0.084 

.  3333 

.3407 

-.0074 

0.012 

.1200 

.1434 

-.0234 

0.222 

1.0000 

.9456 

.0544 

2.763 

.9400 

.9214 

.0187 

0.240 

.9592 

.8786 

.0806 

2.983 

.7083 

.7937 

-.0854 

2.135 

.3750 

.5888 

-.2138 

9.065 

.3265 

.3629 

-.0364 

0.281 

.9600 

.9637 

-.0037 

0.019 

.9348 

.9020 

.0328 

0.559 

.7917 

.7760 

.0157 

0.068 

.5870 

.5663 

.0207 

0.0B0 

.4783 

.2997 

.1786 

6.991 

.1489 

.1403 

.0086 

0.029 

Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

7 

1655.7 

Departure  from  Model 

SSE 

9 

27.125 

TOTAL 

SST 

16 

1682.8 

2 

Coefficient  of  Determination  R  =  .98388 
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3:  “3  =  0  ;  ol  -  3  (Third  and  Fourth  Standardized  Central  Moments)  ; 

=  4.874  ;  p  =  6.158  ;  Five  Parameters  to  be  Estimated 


Parameter  Estimates 
r  =  7.156 
ax  =  0.689 
B1  =  0.286 
a2  =  2.702 
B2  =  0.136 


Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

4 

1673. 8 

Departure  from  Model 

S3E 

12 

27. 130 

TOTAL 

GST 

3  6 

1700.9 

Coefficient  of  Determination  R  =  .98405 
Residual  Analysis 


Pi 

P, 

1 

Residual 

Chi-square 

.8800 

.9124 

-.0324 

0.657 

.8571 

.8020 

.0551 

0.937 

.5217 

.5428 

-.0211 

0.083 

.3333 

.3405 

-.0072 

0.011 

.1200 

.1433 

-.0233 

0.220 

1.0000 

.9461 

.0539 

2.737 

.9400 

.9220 

.0180 

0.225 

.9592 

.8796 

.0796 

2.931 

.7083 

.7951 

-.0868 

2.222 

.3750 

.5909 

-.2159 

9.257 

.3265 

.3649 

-.0384 

0.311 

.9600 

.9641 

-.0041 

0.024 

.9348 

.9029 

.0319 

0.533 

.7917 

.7775 

.0142 

0.056 

.5870 

.5683 

.0187 

0.066 

.4783 

.3014 

.1769 

6.B37 

.1489 

.1412 

.0077 

0.023 

4:  Takahasi  Burr  -  r=  0;ct3  =  0;a4=3  (Third  and  Fourth 
Standardized  Central  Moments);  =  4.874  ; 

p  =  6.158  ;  Four  Parameters  to  be  estimated 


Parameter  Estimates 
ax  =  0.743 
B  =  0.283 
a2  *  2.458 
B2  =  0.145 


Chi-square  Analysis  Table 


source 

djf^ 

Due  to  Model 

SSR 

3 

Departure  from  Model 

SSE 

13 

TOTAL 

SST 

16 

Coefficient  of  Determination  R 


Residual  Analysis 


1483.8 

34.593 

1518.4 

.97722 


P  . 

P. 

Chi-square 

Residual 

1 

.8800 

.9176 

-.0376 

0.935 

.8571 

.8128 

.0443 

0.631 

.5217 

.5617 

-.0400 

0.298 

.3333 

.3600 

-.0267 

0.149 

.1200 

.1571 

-.0371 

0.519 

1.0000 

.9596 

.0404 

2.019 

.9400 

.9391 

.0009 

0.001 

.9592 

.9009 

.0583 

1.867 

.7083 

,8202 

-.1119 

4.074 

.3750 

.6113 

-.2363 

11.284 

.3265 

.3706 

-.0441 

0.408 

.9600 

.9052 

.0548 

1.750 

.9348 

.8417 

.0931 

2.991 

.7917 

.7370 

.0548 

0.742 

.5870 

.5635 

.0235 

0.103 

.4783 

.3025 

.175  8 

6.734 

.1489 

.1343 

.0147 

0.087 

APPENDIX  V 


Data  of  J.  T.  Martin  (1942)  :  The  toxicities  to  Macros iphonie 1 la 
sanborni  of  rotenone,  Sl-elliptone ,  and  of  a  mixture.  Tests  of  7  July  1941. 
Fivefold  replication.  Results  one  day  after  spraying.  Medium  0.5% 
saponin,  containing  5%  alcohol.  Tattersfield  apparatus.  The  data  as 
described  by  Martin,  the  translated  data,  and  the  analyses  of  this  set 
of  data  (data  set  five)  are  in  this  appendix. 
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BIVARj  ATE  LOGIT  ANALYSIS 


Parameter  Estimates 
aQ  =  -1.000 
a1  -  -1.611 
=  -3.509 
a2  =  -3.136 
B2  =  -3.145 


Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

4 

675.39 

Departure  from  Model 

SSE 

7 

30.169 

TOTAL 

SST 

11 

705.56 

2 

Coefficient  of  Determination  R  =  .95724 


Residual  Analysis 


Pi 

P. 

1 

Residual 

Chi-square 

1.0000 

.9240 

.0760 

2.303 

.8500 

.8159 

.0341 

0.  310 

.  3750 

.5164 

-.1414 

3.844 

.0316 

.1510 

-.0694 

1.841 

.9216 

.9221 

-.0005 

0.000 

.9184 

.8275 

.0909 

2.836 

.  3878 

.5095 

-.1217 

2.904 

.0652 

.2133 

-.1481 

6.010 

.9200 

.8549 

.0651 

1.707 

.6667 

.5743 

.0924 

1.678 

.3182 

.2141 

.1041 

2.835 

.1020 

.0441 

.0579 

3.902 

A 

'< 

i 
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BIVARIATE  BURRIT  ANALYSES 

Is  race  -  Vi  nVv^  D2r2?'“t“rc  tC  fee  Estimated 


&l  =  0.224 


B1  = 


0.373 


a2  =  -0.483 
B„  =  0.260 


Farameter 

Estimates 

Chi-square 

Analysis 

Table 

r  = 

7.345 

source 

d.f . 

b,  = 

6.078 

Due  to  Model 

SSR 

7 

1081.9 

*2  = 

D  = 

5.388 

6.345 

Departure  from  Model  SSE 

TOTAL  SST 

4 

11 

29.221 

1111.1 

Coefficient  of  Determination  R 


.97370 


Residual  Analysis 


pi 

P. 

X 

Residual 

Chi-square 

1.0000 

.9691 

.0309 

0.893 

.8500 

.8526 

-.0025 

0.002 

.  3750 

.4718 

-.0967 

1.  803 

.0816 

.0896 

-.0080 

0.039 

.9216 

.9361 

-.0145 

0.180 

.9184 

.8434 

.0750 

2.087 

.  3878 

.5585 

-  1707 

5.794 

.0652 

.  2823 

-.2171 

10.704 

.9200 

.8784 

.0416 

0.811 

.6667 

.5871 

.0796 

1.254 

.3182 

.2258 

.0924 

2.151 

.1020 

.0460 

.0560 

3.503 
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2:  Takahasi  Burr  -  r  =  0  ;  Seven  Parameters  to  be  Estimated 


Parameter  Estimates 
bL  =  9.296 
b2  =  7.696 
p  =  77.198 
ax  =  1.655 
B  =  0.183 
a2  =  0.700 
B2  =  0.142 


f!hi  —  smiarA  inalwoie  pp3.bl's 

— . .  - ..  —  — .  •-  *  r  ? ■=■ . 


source 

d.  f . 

Due  to  Model 

SSR 

6 

1508.3 

Departure  from  Model 

SSE 

5 

31.901 

TOTAL 

SST 

11 

1540.2 

Coefficient  of  Determination 


.97929 


Residual  Analysis 


Pi 

P, 

1 

Residual 

Chi-square 

1.0000 

.9804 

,0196 

0.559 

.8500 

.8623 

-.0123 

0.051 

.3750 

.4921 

-.1171 

2.631 

.0816 

.1317 

-.0501 

1.074 

.9216 

.9509 

-.0293 

0.935 

.9184 

.8434 

.0751 

2.089 

.3878 

.5262 

-.1384 

3.  767 

.0652 

.2614 

-.1962 

9.171 

.9200 

.7874 

.1326 

5.253 

.6667 

.5364 

.1304 

3.280 

.3182 

.2376 

.0806 

1.577 

.1020 

.0602 

.0418 

1.515 
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3:  a ^  =  ;  a 4  =  3  (Third  and  Fourth  Standardized  Central  Moments)  ; 

=  4.874  ;  p  =  6.153  ;  Five  Parameters  to  be  Estimated 

Parameter  Estimates 
r  =  7.158 

=  -0.085 
B  =  0.414 

a2  =  -0.588 
D2  =  0.260 


Residual  Analysis 


fi 

P. 

1 

Residual 

Chi-square 

1.0000 

,9598 

.0402 

1.174 

.8500 

.8426 

.0074 

0.016 

.3750 

.4811 

-.1061 

2.166 

.0816 

.0925 

-.0109 

0.069 

.9216 

.9181 

.0035 

0.009 

.9184 

.8217 

.0697 

3.126 

.3878 

.5500 

-.1622 

5.205 

.0652 

.2887 

-.2235 

11.188 

.9200 

.8762 

.0438 

0.883 

.6667 

.5942 

.0725 

1.046 

.3182 

.2333 

.0849 

1.773 

.1020 

.0504 

.0516 

2.723 

Chi-square  Analysis  Table 


source 

d.  f. 

Due  to  Model 

SSR 

4 

862.11 

Departure  from  Model 

SSE 

7 

29.379 

TOTAL 

SST 

11 

891.49 

2 

Coefficient  of  Determination  R  =  .96705 
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4:  Takahasi  Burr  -  r=  0;ii  =0;a  =3  (Third  and  Fourth 

Standardized  Central  Moments);  =  4.874  ; 

p  =  6.158  ;  Four  Parameters  to  be  Estimated 

Parameter  Estimates 
a  =  -0.080 
B  =  0.418 

a2  =  -0.488 

b  =  0.254 

2  2 

Coefficient  of  Determination  R  =  .96115 


Chi-square  Analysis  Table 


source 

1 

d.f. 

- 1 

Due  to  Model 

SSR 

3 

905.59 

Departure  from  Model 

SSE 

8 

36.603 

TOTAL 

SST 

11 

942.19 

Residual  Analysis 


Pi 

P. 

1 

Residual 

Chi-square 

1.0000 

.9659 

.0341 

0.989 

.8500 

.8588 

-.0088 

0.025 

.  3750 

.5037 

-.1287 

3.179 

.0816 

.099  3 

-.0177 

0.172 

.9216 

.9207 

.0009 

0.001 

.9184 

.8294 

.0890 

2.741 

.3878 

.5697 

-.1819 

6.615 

.0652 

.3118 

-.2466 

13.032 

.9200 

.7923 

.1277 

4.954 

,6667 

.5576 

.1091 

2.317 

.3182 

.2484 

.0698 

1.147 

.  1020 

.0611 

.0409 

1.430 

APPENDIX  VI 


Data  of  J.  T.  Martin  (1942} :  The  toxicities  to  Macros iphonle 11a 
sanborni  of  rotenone,  £-a-toxicarol ,  and  of  a  mixture.  Tests  of 
24  September  1941.  Fivefold  replication.  Results  one  day  after  spraying. 
Medium  0.5%  3aponin,  containing  5%  of  alcohol.  Tattersfield  apparatus. 

The  data  as  described  by  Martin,  the  translated  data,  and  the  analyses  of 
this  set  of  data  (data  set  six)  are  in  this  appendix. 
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DATA  AS  DESCRIBED  BY  MARTIN 


TRANSLATED  DATA 


CONCENTRATIONS  {mg . /I . ) 

Ni 

Number  of 
Insects 
Used 

Percent 

Mortality 

Pi 

X(l) 

Rotenone 

X(2) 

5,-a-Toxicarol 

1.06 

0.0 

51 

100.0 

1.0000 

0.85 

0.0 

43 

97.9 

.9792 

0.64 

0.0 

43 

93.8 

.9375 

0.42 

0.0 

43 

62.5 

.6250 

0.21 

0.0 

43 

12.5 

.1250 

0.0 

9.75 

49 

100.0 

1.0000 

0.0 

7.80 

43 

97.9 

.9792 

0.0 

5.85 

52 

98.1 

.9808 

0.0 

3.90 

49 

87.7 

.8776 

0.0 

1.95 

48 

50.0 

.5000 

0.53 

4.88 

48 

100.0 

1.0000 

0.42 

3.90 

48 

100.0 

1.0000 

0.32 

2.93 

49 

89.8 

.8980 

0.21 

1.95 

50 

82.0 

.8200 

0.11 

0.98 

50 

30.0 

.3000 

<  Q. 
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BIVARIATE  NORMIT  ANALYSIS 

2:  Using  Relative  Change  in  Unweighted  Sum  of  Squares 
Due  to  Error  as  Part  of  the  Convergence  Criteria 


Parameter  Estimates 
a1  «  2.279 

Bx  =  2.063 

a2  =  -0.589 
B2  =  1.290 

=  -0.9999 


Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

4 

.  19" 101 8 

Departure  from  Model 

SSE 

10 

•78*1015 

TOTAL 

SST 

14 

.19‘1018 

2 

Coefficient  of  Determination  R  =  .99592 
Unweighted  Sum  of  Squares  Due  to  Error  =  .03127 


Residual  Analysis 


Pi 

P. 

X 

Residual 

Chi-square 

1.0000 

.9918 

.0082 

0.422 

.9792 

.9740 

.0052 

0.050 

.9375 

.9129 

.0246 

0.367 

.6250 

.6878 

-.0628 

0.882 

.1250 

.1736 

-.0496 

0.790 

1.0000 

.9906 

.0094 

0.465 

.9792 

.9808 

-.0016 

0.006 

.9808 

.9545 

.0263 

0.828 

.8776 

.8784 

-.0008 

0.000 

.5000 

.6074 

-.1074 

2.321 

1.00C 

1.0000 

.0000 

0.000 

1.0000 

1.0000 

.0000 

0.000 

.8980 

1.0000 

-.1020 

0.78* 1015 

.8200 

.7810 

.0390 

0.445 

.3000 

.2806 

.0194 

0.928 
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BIVARIATE  BURRIT  ANALYSES 

i:  General  Case  —  EiyuL  rcuoucteib  Lu  uc  ^alimauQu 

Parameter  Estimates 
r  =  8.817 
b  =  5.045 
b2  =  5.469 
p  =  7.826 
a  =  2.869 
B  =  0.339 
a2  =  3.055 
B2  =  0.184 


Residual  Analysis 


pi 

P. 

1 

Residual 

Chi-square 

1.0000 

.9948 

.0052 

0.266 

.9792 

.9797 

-.0005 

O.OOi 

.9375 

.9144 

.0231 

0.327 

.6250 

.6434 

-.0184 

0.070 

.1250 

.1202 

.0049 

0.011 

1.0000 

.9932 

.0068 

0.335 

.9792 

.9851 

.0059 

0.113 

.9808 

.9612 

.0196 

0.536 

.8776 

.8850 

-.0074 

0.026 

.5000 

.6025 

-.1025 

2.107 

1.0000 

.9994 

.0006 

0.030 

1.0000 

.9920 

.0080 

0.385 

.8980 

.9355 

-.0375 

1.140 

.8200 

.6866 

.1334 

4.17.5 

.  3000 

.2734 

.0266 

0.179 

Chi-square  Analysis  Table 


source 

d.  f . 

Due  to  Model 

SSR 

7 

67305 

Departure  from  Model 

SSE 

7 

9.6604 

TOTAL 

SST 

14 

67314 

2 

Coefficient  of  Determination  R  =  .99986 
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3:  a.  =  0  ;  a.  =  3  (Third  and  Fourth  Standardized  Central  Moments)  ; 
3  4 

=  bj  =  4.874  ;  p  =  6.158  ;  Five  Parameters  to  De  estimated 


Parameter  Estimates 
r  =  5.265 
=  2.907 
B  =  0.357 
a2  =  3.215 
B2  =  0.179 


Residual  Analysis 


pi 

Pi 

Residual 

Chi-square 

1.0000 

.9943 

.0057 

0.294 

.9792 

.9808 

-.0016 

0.006 

.9375 

.9267 

.0108 

0.083 

.6250 

.  6°  ^4 

-.0684 

1.056 

.1250 

.1561 

-.0311 

0.353 

1.0000 

.9826 

.0174 

0.869 

.9792 

.9697 

.0095 

0.148 

.9808 

.9397 

.0411 

1.549 

.8776 

.8637 

.0139 

0.081 

.5000 

.6229 

-.1229 

3.086 

1.0000 

.9980 

.0020 

0.095 

1.0000 

.9868 

.0132 

0.642 

.8980 

.9297 

-.0317 

0.755 

.8200 

.7131 

.  1069 

2.792 

.3000 

.3326 

.0326 

0.240 

Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

4 

23681 

Departure  from  Model 

SSE 

10 

12.047 

TOTAL 

SST 

14 

23693 

2 

Coefficient  of  Determination  R  =  ,99949 
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4:  Takahasi  Burr  -  r  =  0  ;  =  0  ;  =  3  (Third  and  Fourth 


C 4- 


A  W-— \  - 

-S..4W4.W*  WWiMCUWlJ/  / 


“2 


p  =  6.158  i  Four  Parameters  to  be  Estimated 


Parameter  Estimates 
a±  =  2.741 
B  -  0.357 
a2=  3.171 
B2  =  0.189 


Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

3 

6862.7 

Departure  from  Model 

SSE 

11 

23.480 

TOTAL 

SST 

14 

6886.2 

2 

Coefficient  of  Determination  R  =  .99659 


Residual  Analysis 


Pi 

P. 

1 

Residual 

Chi-square 

1.0000 

.9857 

.0143 

0.741 

.9792 

.9567 

.0225 

0.586 

.9375 

.8604 

.0771 

2.374 

.6250 

.5544 

.0706 

0.969 

.1250 

.086  3 

.0387 

0.911 

1.0000 

.9914 

.0086 

0.423 

.9792 

.9840 

-.0048 

0.070 

.9808 

.9648 

.0160 

0.393 

.8776 

.9088 

-.0312 

0.575 

.5000 

.6934 

-.1934 

8.444 

1.0000 

.9773 

.0227 

1.114 

1.0000 

.9479 

.0521 

2.640 

.8980 

.8801 

.0179 

0.148 

.8200 

.7154 

.1046 

2.685 

.3000 

.3814 

-.0814 

1.405 

AFPENDIX  VI I 

Data  of  Ashford  and  Smith  (1964) s  Exposure  to  dust  and  prevalence 
of  pneumoconiosis  for  groups  of  mine  workers.  The  data  as  described  by 
Ashford  and  Smith,  the  computed  p^  ,  and  the  analyses  on  this  set  of  data 
(data  set  seven)  are  in  this  appendix. 


<  a 
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BIVARIATE  NUKM1T  ANALYSIS 


Parameter  Estimates 
a  =  -2.818 
B  =  0.937 

a2  =  -2.446 
B  =  0.501 

2  2 

=  -0.320  Coefficient  of  Determination  R  =  .85095 


Chi-aquare  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

4 

217.75 

Departure  from  Model 

SSE 

35 

38.141 

TOTAL 

SST 

39 

255.89 

Residual  Analysis 


Pi 

P. 

1 

Residual 

Chi-square 

.0222 

.0195 

.0027 

0.053 

.1111 

.0800 

.0311 

0.236 

.0625 

.1236 

-.0611 

0.551 

.1765 

.1579 

.0186 

0.044 

.1429 

.1834 

-.0405 

0.153 

.2500 

.2200 

.0300 

0.063 

.2273 

.2451 

-.0179 

0.038 

.2258 

.2703 

-.0445 

0.311 

.2000 

.3030 

-.1030 

1.256 

.2941 

.3091 

-.0150 

0.018 

.0875 

.1488 

-.0613 

2.  375 

.1000 

.2025 

-.1025 

0.651 

.3571 

.2719 

.0852 

0.514 

.2500 

.2886 

-.0386 

0.058 

.5238 

.3120 

.2118 

4.390 
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Residual  Analysis — Continued 


pi 

P. 

1 

Residual 

Chi -square 

.3571 

.3564 

.0007 

0.000 

.5385 

.  3674 

.1710 

1.636 

.7000 

.4202 

.2796 

3.  209 

.2676 

.  3106 

-.0429 

0.612 

.1250 

.3586 

-.2336 

1.898 

.5000 

.3897 

.1103 

0.205 

.2857 

.4020 

-.1163 

0.394 

.5000 

.4515 

.0485 

0.057 

.6000 

.5213 

.0787 

0.248 

.5000 

.5533 

-.0533 

0.046 

.5000 

.4362 

.0638 

1.753 

.4000 

.4631 

-.0631 

0.080 

.4000 

.5200 

-.1200 

0.288 

.6667 

.5945 

.0722 

1.298 

.2500 

.5748 

-.3248 

1.727 

.5862 

.5361 

.0501 

0.586 

0.0000 

.6290 

-.6290 

5.086 

.2500 

.6632 

-.4132 

3.058 

.  3333 

.6666 

-.3333 

1.500 

.6515 

.6178 

.0338 

0.318 

.7500 

.7145 

.0355 

0.025 

1.0000 

.7447 

.2553 

0.686 

.6667 

.6798 

-.0132 

0.026 

.5500 

.7302 

-.1602 

3.295 

.8000 

.7746 

.0254 

0.037 
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BIVARIATE  BURRIT  ANALYSES 

li  General  Cape  -  Eight  Parameters  to  be  Estimated 


Parameter  Estimates 


Chi-square  Analysis  Table 


source 

m 

Due  to  Model 

SSR 

7 

217.30 

Departure  from  Model 

SSE 

32 

38.293 

TOT  AIi 

SST 

39 

255.59 

0.147 

3.533 

0.082 


Coefficient  of  Determination  R  =  .85018 


Residual  Analysis 


p. 

1 

1  P. 

l 

Residual 

Chi-square 

0222 

.0183 

.0039 

0.113 

1111 

.0734 

.0377 

0.377 

0625 

.1170 

-.0545 

0.459 

1765 

.1528 

.0236 

0.073 

1429 

.1809 

-.0381 

0.137 

2500 

.2208 

.0292 

0.060 

2273 

.2489 

-.0217 

0.055 

2258 

.2774 

-.0516 

0.412 

2000 

.3145 

-.1]45 

1.520 

2941 

.3223 

-.0282 

0.062 

0875 

.1500 

-.0625 

2.448 

1000 

.1997 

-.0997 

0.622 

3571 

.2711 

.0860 

0.524 

2500 

.2908 

-.0408 

0.064 

5238 

.3171 

.2067 

4.142 

Residual  Analysis — Contvnued 


1 


2:  TaXahasi  Burr  -  r  =  0  ;  Seven  Parameters  to  be  Estimated 


Parameter  Estimates 


Chi-square  Analysis  Table 


Residual  Analysis — Continued 


- r 

Pi 

- r 

P. 

i 

1 

Residual 

- 1 

Chi-square 

.3571 

.  3507 

.0064 

0.003 

.5385 

.  3554 

.1831 

1.903 

.7000 

.3960 

.3040 

3.865 

.2676 

.3165 

-.0489 

0.784 

.1250 

.  3683 

-.2433 

2.036 

.5000 

.3937 

.1064 

0.190 

.2857 

.4001 

-.1144 

0.  382 

.5000 

.4364 

.06  36 

0.099 

.6000 

.4876 

.1125 

0.506 

.5000 

.  5059 

-.0059 

0.001 

.5000 

.4342 

.0658 

1.867 

.4000 

.4615 

-.0615 

0.076 

.4000 

.5057 

-.1057 

0.223 

.6667 

.5663 

.1004 

0.246 

.2500 

.5428 

-.2928 

1.382 

.5862 

.5276 

.0586 

0.798 

0.0000 

.6110 

-.6110 

4.712 

.2500 

.6388 

-.3888 

2.620 

.3333 

.6218 

-.2885 

1.062 

.6515 

.6056 

.0460 

0.584 

.7500 

.6760 

.0740 

0.100 

1.0000 

.7149 

.2851 

0.798 

.6667 

.6666 

.0001 

0.000 

.5500 

.7176 

-.1676 

2.773 

.8000 

.76  32 

.0368 

0.075 
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3:  a3  =  0  >  Cl4  =  3  (Third  and  Fourth  Standardized  Central  Moments)  ; 
b^  =  =  4.874  ;  p  =  6.158  ;  Five  Parameters  to  be  Estimated 


Parameter  Estimates 
r  =  4.742 
ax  =  1.259 
B  -  0.151 
a2  =  3.405 
B2  =  0.078 


Residual  Analysis 


pi 

P. 

1 

Residual 

Chi-square 

.0222 

.0211 

.0011 

0.008 

.1111 

.0905 

.0206 

0.093 

.0625 

.1340 

-.0715 

0.704 

.1765 

.1670 

.0094 

0.011 

.1429 

.1913 

-.0485 

0.212 

.2500 

.2261 

.0239 

0.039 

.2273 

.2499 

-.0226 

0.060 

.2258 

.2740 

-.0481 

0.361 

.2000 

.3061 

-.1061 

1.326 

.2941 

.3103 

-.0162 

0.021 

.0875 

.1528 

-.0653 

2.634 

.1000 

.2139 

-.1139 

0.772 

.3571 

.2820 

.0752 

0.391 

.2500 

.2978 

-.0478 

0.087 

.5238 

.3199 

.2039 

4.013 

Chi-square  Analysis  Table 


source 

d.f. 

Due  to  Model 

SSR 

4 

218.53 

Departure  from  Model 

SSE 

35 

38.338 

TOTAL 

SST 

39 

256.87 

2 

Coefficient  of  Determination  R  =  .85075 


Residual  Analysis — Continued 


pi 

"i 

Residual 

Chi-squarc 

.3571 

.3631 

-.0060 

0.002 

.5  385 

.3660 

.1725 

1.666 

.7000 

.4127 

.2873 

3.405 

.2676 

.3187 

-.0511 

0.854 

.1250 

.3793 

-.2543 

2.198 

.5000 

.4103 

.0897 

0.133 

.2857 

.4198 

-.1341 

0.517 

.5000 

.4615 

.0385 

0.036 

.6000 

.5201 

.0800 

0.256 

.5000 

.5420 

-.0419 

0.028 

.5000 

.4353 

.0647 

1.806 

.4000 

.4749 

-.0749 

0.112 

.4000 

.5294 

-.1294 

0.  336 

.6667 

.5964 

.0703 

0.123 

.2500 

.5760 

-.3260 

1.740 

.5862 

.5273 

.0589 

0.807 

0.0000 

.6278 

-.6278 

5.059 

.2500 

.6612 

-.4112 

3.018 

.3333 

.6598 

-.3265 

1.425 

.6515 

.6039 

.0477 

0.626 

.7500 

.7086 

.0414 

0.033 

1.0000 

.7396 

.2604 

0.704 

.6667 

.6638 

.0029 

0.001 

.5500 

.7139 

-.1639 

2.632 

.8000 

.7605 

.0395 

.  . .  - . . 

0.086 
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4:  Takahasi  Burr  -r=0;a3=0;a4=3  (Third  and  Fourth 
Standardized  Central  Moments) ;  b^  =  =  4.874  ; 

p  =  6.158  ;  Four  Parameters  to  be  Estimated 
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Residual  Analysis — Cnnti.nueri 


pi 

P 

i 

Residual 

Chi-square 

.3571 

.3542 

.0029 

0.001 

.5385 

.  3581 

.1804 

1.841 

.7000 

.3983 

.3017 

3.799 

.2676 

.3210 

-.0534 

0.930 

.1250 

.3735 

-.2485 

2.110 

.5000 

.3984 

.1016 

0.172 

.2857 

.4045 

-.1188 

0.410 

.5000 

.4399 

.0601 

0.088 

.6000 

.4899 

.1101 

0.485 

.5000 

.5075 

-.0075 

0.001 

.5000 

.4363 

.0637 

1.750 

.4000 

.4645 

-.0645 

0.084 

.4000 

.5079 

-.1079 

0.233 

.6667 

.5669 

.0997 

0.243 

.2500 

.5441 

-.2941 

1.395 

.5862 

.5271 

.0591 

0.814 

0.0000 

.6096 

-.6096 

4.685 

.2500 

.6369 

-.3869 

2.588 

.3333 

.6209 

-.2876 

1.054 

.6515 

.6027 

.0489 

0.658 

.7500 

.6734 

.0766 

0.107 

1.0000 

.7108 

.2892 

0.814 

.  6667 

.6619 

.0048 

0.003 

.5500 

.7115 

-.1615 

2.542 

.8000 

.  7563 

.0437 

0.104 

APPENDIX  VIII 


Listings  of  FORTRAN  subroutines  used  .in  evaluating  the  functions, 
partial  derivatives ,  and  the  weights  are  in  this  appendix. 
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o  o 


1. 


Bivariate  Normal  Subroutines 


DOUBLE  PRECISION  FUNCTION  F(X,A) 

C  BIVARIATE  NORMAL  FUNCTION 
DIMENSION  X(1),A(1) 

DOUBLE  PRECISION  X , A,AA,B,R, S ,GOFU,THA  ,C,D 
AA=A (1) +A (2) *X  (1) 

B=A (3)+A(4)*X(2) 

IF  (DABS (A (5) ) .LE. 0. 9999D+00)  GO  TO  16 
A(5)«DSIGN(0.9999D+00,A(5) ) 

16  R=A (5) 

S^DSQRT (1 . -R*R) 

C=GOFU (AA) 

D=GOFU(B) 

FsC+D-THA (AA,B/AA) -THA(B,AA/B) 

1+THA (AA, (B-R*AA)/(AA*S) )+THA(B, (AA-R*B) / (B*S) ) -C*D 
RETURN 
END 


SUBROUTINE  PD(X,A,FXA,P) 

DIMENSION  X(l) ,A(1) ,P (1) 

DOUBLE  PRECISION  X,A,AA,R,S,B,FXA,P ,WATE,GOFU,GPRIME,C,D 
A1=ALPHA1 , A2=BETA1 , A3=ALPHA2 , A4=BETA2 , A5=RHO 
BIVARIATE  NORMAL  PARTI ALS 
AA=A  (1)'+A  (2)  *X  (1) 

B=A(3) +A (4) *X (2) 

R-A(5) 

S=DSQRT (l.-R*R) 

C=GPRIME(AA) 

D=  (B-R*AA) /S 
P  (.l)=C*  (l.-GOFU(D)  ) 

P(2)-X(1)*P(1) 

P (3) =GPRIME (B) * (1. -GOFU ( (AA~R*B) /S) ) 

P  (4)  =X  (2)  *P  (3) 

P (5 ) =-C*GPRIME (D) /S 

RETURN 

END 
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DOUBLE  PRECISION  FUNCTION  GPRIME(X) 
DOUBLE  PRECISION  X,A1,A2,A3 
A1  =X*XM-0.5D  0) 

IF  (A1.LE.-150)  GO  TO  15 
A2=DEXP (Al) 

A3  =  . 398942280401433D+00 
GPRIME  =  A2*A3 
GO  TO  16 

15  GPRIME-0 

16  RETURN 
END 


DOUBLE  PRECISION  FUNCTION  GOFU (U) 

DIMENSION  Y (160) 

DOUBLE  PRECISION  U ,GU ,X ,XSET, DELTA, GPR,DELK ,DELTAX ,Y ,SUM, TOP 
IF (Y  (17) - .6914624612 7D+00) 2,5,2 


Y< 

1) 

= 

.  50000000000D+00 

Y< 

2) 

= 

. 39894228040D+00 

Y< 

3) 

= 

.OOOOOOOOOOOD+OO 

Y  ( 

4) 

= 

-.66490380066D-01 

Y< 

5) 

= 

.OOOOOOOOOOOD+OO 

Y( 

6) 

= 

. 99735570100 D-02 

Y  ( 

7) 

= 

.OOOOOOOOOOOD+OO 

Y  ( 

8) 

= 

- . 11873282155D-02 

Y  ( 

9) 

- 

. 59870632568D+00 

Y  ( 

10) 

= 

.  38666811680 D+OO 

Y( 

11) 

= 

- . 48333514600D-01 

Y  ( 

12) 

= 

-. 604 1689 3250D-01 

Y  ( 

13) 

= 

. 11831641595D-01 

Y  ( 

14) 

= 

. 84709519078D-02 

Y  ( 

15) 

= 

- . 19305085421D-02 

Y  ( 

16) 

= 

-. 93949992204D-03 

Y  ( 

17) 

= 

.  69146246127D+00 

Y  ( 

18) 

= 

. 352065 32676D+00 

Y  ( 

19) 

- . 88016  33 1691D-01 

Y  ( 

20) 

= 

-. 44008165845D-01 

Y  ( 

21) 

= 

. 20170409 346D-01 

Y( 

22) 

= 

. 4584183942 3D-02 

Y  ( 

23) 

= 

- . 30714032413D-02 

Y  ( 

24) 

= 

- .  32635023779D-03 

Y  ( 

25) 

- 

. 7733726476 3D+00 

Y  ( 

26) 

= 

. 30113743216D+00 

Y  ( 

27) 

= 

- .  11292653706D+00 

Y  ( 

28) 

- 

2 19 5 79 37 76 ID- 01 

Y  ( 

29) 

= 

. 22938202840D-01 

Y( 

30) 

= 

14703976180D-03 
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Y< 

31) 

= 

-.30400470751D-02 

Y( 

32) 

. 34322406 302D-03 

Y( 

33) 

= 

.841344746060+00 

Y  ( 

34) 

- 

•24197072452D+00 

Y  ( 

3b) 

3 

-.120995362260+00 

Y  ( 

36) 

= 

.000000000000+00 

Y  ( 

37) 

= 

.  20164227043D-01 

Y  ( 

38) 

= 

- . 40328454087D-02 

Y( 

39) 

s* 

- . 20 16422 704 3D-02 

Y( 

40) 

= 

. 7681610 3022D-03 

Y( 

41) 

= 

. 89435022633D+00 

Y( 

42) 

= 

.182649085 39D+00 

Y  ( 

43) 

- 

11415567837D+00 

Y< 

44) 

= 

.171233517550-01 

Y( 

45) 

B 

.136748989 7 lD-01 

Y  ( 

46) 

a 

- .  598722 75061D-02 

Y( 

47) 

= 

-.575980799060-03 

Y( 

48) 

= 

.815618893410-03 

Y( 

49) 

= 

. 93319279874D+00 

Y( 

50) 

= 

. 12951759567D+00 

Y  ( 

51) 

= 

-. 97138196750D-01 

Y( 

52) 

= 

. 26982832430D-01 

Y  ( 

53) 

= 

.60711372 969D-02 

Y( 

54) 

= 

-. 586876605 36D-02 

Y( 

55) 

a 

.  65770654049D-03 

Y( 

56) 

= 

.557725509610-03 

Y( 

57) 

= 

« 95994084 314D+0Q 

Y( 

58) 

= 

. 8627731882 7D-01 

Y( 

59) 

= 

- .  7549265 3974D-01 

Y( 

60) 

= 

. 29657828346D-01 

Y( 

61) 

= 

-.393190906110-03 

Y( 

62) 

= 

-.431105743480-02 

Y  ( 

63) 

a 

. 13098172060D-02 

Y  ( 

64) 

= 

.185766821700-03 

Y  ( 

65) 

= 

. 97724986 805D+00 

Y  ( 

66) 

= 

.5 39909665 13D-01 

Y  ( 

67) 

= 

-.539909665130-01 

Y( 

68) 

= 

.  269954 8325 7D-01 

Y( 

69) 

= 

- . 4 499 2 4 72 0 9 4D- 02 

y  ( 

70) 

= 

-.224962360470-02 

Y( 

71) 

= 

. 13497741628D-02 

Y  ( 

72) 

- 

- . 1178 3 74269 lD-0  3 

73) 

= 

.987775527350+00 

Y  ( 

74) 

= 

. 31739651836D-01 

Y  ( 

75) 

= 

-.  35707108315D-01 

Y  ( 

76) 

= 

.  21490389264D-01 

Y( 

77) 

= 

-.61371592417D-02 

Y( 

78) 

= 

- .  46183673081D-03 

Y( 

79) 

= 

.99147667294D-03 

Y( 

80) 

= 

-.263708367400-03 

Y(  81)  = 

Y  (  82)  = 

Y  (  83)  = 
Y(  84)  = 
Y<  85)  = 
Y(  86)  = 
Y (  87)  = 
Y(  88)  = 
Y<  89)  = 
Y(  90)  = 

Y  {  91)  = 

Y  <  92)  = 

Y  (  93)  = 
Y(  94)  = 
Y (  95)  = 
Y (  96)  = 

Y  (  97)  - 
Y(  98)  = 
Y(  99)  = 

Y  (100)  - 
Y (101)  = 
Y (102)  = 
Y (103)  = 
Y (104)  = 

Y  (105)  = 
Y (106)  = 
Y (107)  - 
Y (108)  = 
Y (109)  « 
Y (110)  = 
Y(lll)  = 

Y  (112 )  = 
Y(113)  = 
Y (114)  = 
Y (115 )  = 
Y (116)  = 

Y  (117)  = 

Y  (118)  = 

Y  (119)  = 
Y (120)  = 

Y  (121)  - 

Y  (122)  = 

Y  (123)  = 

Y  ( 124)  = 
Y (125 )  = 
Y (126 )  = 

Y  (127)  = 

Y  (128)  = 

Y  (129)  = 

Y  (130)  = 


.99379033467D+00 
-17528300494D-01 
2 191 0  37661  7n-m 
. 15337262932D-01 
- . 59340600629D-02 
•66644059169D-03 
•51352442853D-03 
-.26273974729D-03 
.99702023677D+00 
.90935625017D-02 
12503648440D-01 
.99 4608 3986 ID-02 
- . 47539913338D-02 
.  11227826 357D-02 
. 11925680315D-03 
18051548644D-03 
.99 8650 1019 7D+00 
. 44318484120D-02 
66 477726 180D-02 
.  5909] 312160D-02 
-. 332388630 90D-02 
. 11079621030D-02 
-. 11079621030D-03 
- . 84416160226D-04 
.  99942297496D+00 
.  20290480573D-02 
- . 32972030931D-02 
. 32337953413D-02 
20779248659D-02 
. 86558186169D-03 
-. 19180019295D-03 
13995370 141D-04 
. 99976737091D+00 
. 87268269505D-03 
-.1527194716 3D-  02 
. 16362800532D-02 
- . 117721259 38D-0 2 
. 57860680771D-03 
-. 18055895865D-03 
. 21397716502D-04 
. 9999 11582 7 1D+00 
. 352595682 37D-0 3 
-. 661116 90444D-0 3 
. 76763018349D-03 
- . 60946714628D-03 
, 34195583218D-03 
13246010895D-03 
. 30251745009D-04 
. 99996832876D+00 
. 13383022576D-03 
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Y  (131) 
Y (132) 
Y (133) 
Y (134) 
Y (135 ) 

Y  (136) 
Y (137) 
Y (138) 
Y (139) 
Y (140) 
Y ( 141 ) 
Y (142 ) 
Y ( 143) 
Y (144) 
Y (145 ) 
Y (146) 
Y (147) 
Y (148) 
Y (149) 
Y (150 ) 
Y (151) 
Y(152) 
Y (153) 
Y (154) 
Y (155 ) 
Y (156) 
Y (157) 
Y (158) 
Y (159 ) 
Y (160) 


26766045153D-03 
. 33457556441D-03 
289965489 16D-03 
. 18178605667D-03 
- . 82528639221D-04 
. 255 18025 191D-04 
. 99998931147D+00 
.4 77186365 40D-04 
-.10140210265D-03 
. 13569987266D-03 
- . 12728076426D-03 
.8783366 87 23D- 04 
-.45244746779D-04 
. 17013635696D-04 
.999996602330+00 
. 15983741107D-04 
- . 35963417490D-04 
. 51281169385 D-04 
- . 51697412642D-04 
.  38835495972D-04 
-.22233S33626D-04 
. 96697768581D-05 
. 99999698292D+00 
. 50295072886D-05 
-. 11945079810D-04 
.180747918180-04 
-.194729686490-04 
. 15788101444D-04 
-.990251782340-05 
. 48400297796D-05 


5  IF  (U)  11,10,12 

10  GU-.5D+00 
GO  TO  100 

11  X=DABS(U) 

GO  TO  13 

12  X=U 

13  IF (X-7.0D+00) 15,14,14 

14  IF (U) 141, 10, 142 

141  GU=0. 00+00 
GO  TO  100 

142  GU-1. 00+00 
100  GOFU=GU 

RETURN 

15  IF (X-4.87499D+00) 16 ,16,40 

16  XSET=X*4 . 00+00 
XSST=XSET+. 50+00 
I=IFIX  (XSST) 
XSET-DFLOAT(I) 

DELTA=X- (XSET* . 25D+00) 

201  K=I*8+1 
I=K+7 

SUM=0. 0D+00 
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DO  20  J=K, I , 1 
L=K+I-J 
SUM=SUM*DELTA 
SUM=Y (L) +SUM 

20  CONTINUE 
IF(tJ)  21,10,22 

21  GU=1.0D+00  -SUM 
GO  TO  100 

22  GU-SUM 

GO  TO  100 

40  XSET=-X*X 

IF  (XSET.LE.-300)  GO  TO  101 

GPR= (DEXP (XSET* . 5D+0Q ) ) *. 3989 4228040 14 33D+00 

GO  TO  102 

101  GPR=0 

102  DELTA® 1 . 0D+00/X 
SUM® DELTA 
T0P®1 . OD+OO 

41  DELK®TOP/XSET 
DELTAX=DELTA*DELK 

IF  (DABS (DELTA) -DABS (DELTAX) ) 45 ,45 ,43 
43  DELTA=DELTAX 
SUM=SUM+DELTA 

IF (DABS (GPR*DELTA) - . 5D-9) 45,45,42 

42  TOP=TOP+2 . OD+OO 
GO  TO  41 

45  SUM=GPR*SUM 
IF (U)  22,10,21 
END 


DOUBLE  PRECISION  FUNCTION  THA  (HX,AX) 

DOUBLE  PRECISION  HX ,AX , AA,U ,H , A , SUM,C ,DA ,TA ,TX, X, Y , Z  ,GOFU 
DIMENSION  AA  (9)  ,U(9) 

DATA  AA (1) ,AA (9) ,AA(2) ,AA(8) ,AA(3) ,AA(7) ,AA(4) ,AA(6) ,AA (5) /2* . 4063 
17194181E-1 ,2* . 903240^0 347E-1, 2*. 1303053482 , 2* . 15617353852, . 1651196 
2775/, U(l) , U (2 ) ,U (3) ,U(4) ,U(5) ,U(6) ,U(7) ,U(8) ,U(9)/.15919880246E-1, 
3. 81984446 337E-1,.  19331428365,. 3378732883,. 5,  .6621267117 ,. 806685716 
435, .91801555366, .98408011975/ 

H=HX 

A=AX 

IF (DABS (H) .LE.5.77)GO  TO  10 

11  THA=0 . 

RETURN 

10  IF  (DABS  (A)  .  LE .  1 . )  GO  TO  13 

12  H=A*H 

IF  (DABS  (H)  .  LE.  5 . 77)  GO  TO  15 
GO  TO  16 
15  A=l./A 


13  SUM=0. 

DO  61  M=1 ,9 , 1 
DA=1.+A**2*U(M) **2 
C=-.5*H**2*DA 
TA=DEXF  (C) /DA  *  AA (M) 
SUM=SUM+TA 
61  CONTINUE 

TX=A/6.2831853072*SUM 
GO  TO  17 

16  TX-O. 

17  IF  (DABS  (AX) .LE.l. )GO  TO  20 

14  X=GOFU (HX) 

Y=GOFU<H) 

Z=X*Y 

TX=.5*X+.5*Y-Z-TX 

18  IF (AX) 21 ,20 ,20 
21  TX=TX- . 5D+00 
20  THA=TX 

RETURN 

END 
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2.  Bivariate  Logistic  Subroutines 


DOUBLE  PRECISION  FUNCTION  F(X,A) 

C  BIVARIATE  LOGISTIC  (GUMBEL)  FUNCTION 
DIMENSION  X(l) , A (1 ) 

DOUBLE  PRECISION  X ,A ,C,CC , B , BB,D,E , BC ,DE 

IF  (DABS (A (1) ) . GE . 1. OD+OO)  A (1) =DSIGN (1. OD+OO ,A (1) ) 

C=A  (3)  *  (X  (1)  +  A  (2)  ) 

IF  (DABS (C) .GT. 150. OD+OO)  C=DSIGN (150. OD+OO ,C) 
CC=DEXP (C) 

B=A (5) * (X (2)+A (4) ) 

IF  (DABS (B) .GT. 150. OD+OO)  B=DSIGN (150. OD+OO ,B) 
BB=DEXP (B) 

D=l./ (l.+CC) 

E=l. / (l.+BB) 

BC=BB*CC 

DE=D*E 

F=D+E-DE* < 1. +A (1) *BC*DE) 

RETURN 

END 


SUBROUTINE  PD (X ,A, FXA,P ) 

BIVARIATE  LOGISTIC  (GUMBEL)  PARTIALS 

A ( 1 ) = ALP HA  0  ,  A  ( 2 ) =ALPHA1 ,A  (3) =BETA1 , A ( 4 ) =ALPHA2 ,A (5) =BETA2 
DIMENSION  X(l) ,A(1) ,P(1) 

DOUBLE  PRECISION  X ,A ,P ,FXA , WATE,C ,CC ,B ,BB,D ,DD ,E ,EE ,BC ,DE , R,S ,T ,Z 
C=A  (3)  *  (X  (1)  +A  (2)  ) 

IF  (DABS (C) .GT. 150. OD+OO)  C=DSIGN  (150. OD+OO ,C) 

CC=DEXP (C) 

B=A (5)  *  (X  (2)  +A  ( 4)  ) 

IF  (DABS (B).GT. 150. OD+OO!  B=DSIGN (150. OD+OO ,B) 

BB=DEXP (B) 

D=l. /(l.+CC) 

DD=D*D 
E=l. /(l.+BB) 

EE=E*E 

BC=BB*CC 

DE=D*E 

Z=1 . +A (1) *BC*DE 
P  ( 1 ) =-DD*EE*BC 
T=A (1) *P(1) 

R=DD*CC* (E*Z-1. )+T* (l.-D+CC) 

P  (2)  =A  (3)  *R 
P(3)=(X(1)+A(2)  )*R 
S=EE*BB*(D*Z-1, )+T* (l.-E+BB) 

P(4)=A(5)*S 

P  (5)  =  (X  (2 )  +A  (4)  )  *S 

RETURN 

END 
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3.  Bivariate  Burr  Subroutines 


DOUBLE  PRECISION  FUNCTION  F(X,A) 

C  BIVARIATE  BURR  R  (FD)  FUNCTION 

C  A1=R,A2=B1 , A3=B2 , A4=F , A5=ALPHA1, A6=BETA1 , A7=ALPHA2 , A8=BETA2 

DOUBLE  PRECISION  X,A ,B ,BB ,C ,CC ,D,DD 
DIMENSION  X(l) ( A (1) 

IF (A (1) . GT.A(4)+1.0D+00)  A (1)=A (4) +1 . 0D+00 
IF(A(1) .LE.0.0D+00)  A(1)=0.0 
B=A (6) * (X ( 1) +A (5 ) ) 

IF  (B.GT.0.0D+00)  GO  TO  17 
BB=0 

GO  TO  18 

17  IF  (A(2) *DLOG(B) .LT. 150.0)  GO  TO  19 
BB=DEXP (150.0D+00) 

GO  TO  18 

19  IF  (A (2) ) 13 , 14 , 15 

13  BB= 1 . / (B**DABS (A (2) ) ) 

GO  TO  18 

14  BB=1. 0 

GO  TO  18 

15  BB=B**A (2 ) 

18  C=A  (8)  *  (X  (2 )  +A  (7) ) 

IF  (C.GT.0.0D+00)  GO  TO  16 

CC=0 

GO  TO  10 

16  IF (A ( 3 ) *DLOG (C ) . LT .150.0)  GO  TO  20 
CC=DEXP (150 . OD+OO) 

GO  TO  10 

20  IF  (A (3) ) 9 , 11 , 12 

9  CC=1./ (C**DABS (A (3) ) ) 

GO  TO  10 

11  CC=1.0 
GO  TO  10 

12  CC=C**A (3) 

10  D= (1. 0+BB+CC+A (1) *BB*CC) 

IF  (A(4) *DLOG(D) .LT. 150.0)  GO  TO  21 
DD=DEXP ( -150 . OD+OO ) 

GO  TO  22 

21  DD=1 . 0/ (D**A (4) ) 

22  F=l. 0D+00-DD 
RETURN 

END 
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SUBROUTINE  PD (X ,A ,FXA,P) 

BIVARIATE  BURR  R  PARTTALS 

A1=R,  A2=Bl ..  A3-=B2  ,  A4=P  , A5-ALPHA1 ,  A6=BETA1 ,  A7=ALPHA2  ,  A8=BETA2 
DIMENSION  Xfl) ,A(1) ,P(1) 

DOUBLE  PRECISION  X , A ,FXA,P , WATE , B , BB,C ,CC ,D, DD ,E ,EE ,G ,GG,H,HH , I 
IF (A (1) . GT.A (4)+1.0D+00)  A(1)=A(4)+1.0D+00 
IF(A(1) .LE.0.0D+00)  A(1)=0.0 
B=A (6) * (X (1) +A (5) ) 

IF  (B.GT.0.0D+00)  GO  TO  17 
BB=0 

GO  TO  18 

17  IF  (A(2)*DLOG(B) .LT. 150.0)  GO  TO  30 
BB=DEXP(150.0D+00) 

GO  TO  18 

30  IF  (A(2 )) 13, 14,15 

13  BB=1./(B**DABS<A<2) ) ) 

GO  TO  18 

14  BB=1 . 0 
GO  TO  18 

15  BB=B**A (2) 

18  C=A  (8)  *  (X  (2)  +A  (7)  ) 

IF  (C.GT.O.OD+OO)  GO  TO  16 
CC=0 

GO  TO  10 

16  IF (A (3) *DLOG (C) . LT. 150.0)  GOTO  31 
CC=DEXP(150.0D+00) 

GO  TO  10 

31  IF  (A (3) ) 9 , 11 , 12 

9  CC=1. / (C**DABS (A ( 3) ) ) 

GO  TO  10 

11  CC=1.0 
GO  TO  10 

12  CC=C**A (3) 

10  D= (1 . 0+BB+CC+A (1) *BB*CC) 

IF( (A(4)+1.0) *DLOG(D) .LT. 150.0)  GOTO  32 
DD=A (4) *DEXP (-150.0D+00) 

GO  TO  33 

32  DD=A  (4)  *  (  (1.0/D)  **  (A  (4)  +1 . 0)  ) 

33  E=DD*BB 

G=l. 0+A (1) *BB 
GG=1. 0+A (1) *CC 
H=DD*A (2) *GG 
HH=DD*A (3) *G 
P (1) =E*CC 

IF  (B.GT.O.OD+OO)  GO  TO  1 
P  (2)  =0 
GO  TO  2 

1  P  (2) =E*DLOG (B) *GG 

2  IF  (C.GT.O.OD+OO)  GO  TO  3 
P(3)=0 

GO  TO  4 

3  P<3)=DD*CC*DLOG(C)*G 

4  IF  (A(4) *DLOG (D) .LT. 150.0)  GO  TO  34 
P (4)= (DLOG (D) ) * (DEXP (-150.0D+00) ) 
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GO  TO  35 

34  P(4)=DLOG(D)/(D**A(4)) 

35  ir  (B.GT.G.GD-rGG)  GO  TO  23 
P (5) =0 .0 

P (6) =0.0 

GO  TO  24 

23  IP  ((A(2)-1.0)*DLOG(B) .LT. 150.0)  GO  TO  36 
EE=H*DEXP (150. OD+OO) 

GO  TO  8 

36  IF  (A (2) -1.0) 5,6,7 

5  EE=H/ (B**DABS (A (2) -1.0) ) 

GO  TO  8 

6  EE=H 

7  EE=H* (B** (A (2) -1.0)) 

8  P(5)=EE*A  (6) 

P  (6) =EE*X  (1) 

24  IF  (C.GT.O. OD+OO)  GO  TO  25 
P(7)=0.0 

P ( 8 ) =0 . 0 

GO  TO  26 

25  IF  ( (A(3)-1.0) *DLOG (C) .LT. 150,0)  GO  TO  37 
I=HH*DEXP (150. OD+OO) 

GO  TO  22 

37  IF  (A(3)-l. 0)19,20, 21 

19  I=HH/ (C**DABS (A  (3) -1 . 0) ) 

GO  TO  22 

20  I=HH 

GO  TO  22 

21  I=HH* (C** (A(3)-1.0) ) 

22  P  (7)  =I*A  (8) 

P  (8)=I*X (2) 

26  CONTINUE 
RETURN 
END 
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4.  Weight  Subroutine 


DOUBLEPREC I S IONFUNCTIONWATE ( X , FXA ) 
DOUBLEPRECIS I ONX ,FXA 
DIMENSIONX (1) 

WATE=DABS ( X ( 3 ) / ( (1 . 0D+00-FXA) *FXA) ) 

RETURN 

END 
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3  ABSTRACT 

>  The  problem  of  a  mixture  of  two  stimulants  in  a  biological  quantal  assay  is 
investigated  from  a  mathematical  standpoint.  The  basic  assumption  is  made  that 
the  response  region  does  not  depend  on  biological  considerations  -  i.e.,  given  a 
specified  mixture  of  stimulants  z  ,  the  response  region  is  defined  by  the  point 
z'  in  the  p-variate  space  where  there  are  p  stimulants  under  consideration;  instead, 
the  probability  functions,  themselves,  may  take  on  different  forms.  A  general 
form  is  proposed  and  investigated.  Three  analytic  models  (one  utilizing  the  bi¬ 
variate  normal  distribution,  one  a  bivariate  logistic  distribution  developed  by 
Gumbel  (1961)  ,  and  one  a  bivariate  Burr  distribution  developed  by  this  author)  art 
employed  in  this  investigation.  The  investigation  includes  the  analysis  of  data, 
under  the  three  analytic  models,  which  had  been  classified  by  previous  investigators 
as  examples  of  synergistic  action,  simple  similar  action,  independent  action,  and 
additive  accion.  The  residual  analyses  are  included  as  well  as  the  FORTRAN  IV  sub¬ 
routines  used  in  evaluating  the  functions,  the  partial  derivatives  and  the  weights,  i 
The  investigation  lends  some  support  to  the  assumption  of  a  constant  response  ^ 
region  for  a  diversity  of  mixtures  of  stimulants.  The  analytic  model  incorporating 
the  bivariate  Burr  distribution  is  recommended  for  all  cases  unless  the  number  of 
parameters  to  be  estimated  is  a  primary  concern,  in  which  case  the  analytic  model 
utilizing  the  bivariate  normal  distribution  is  recommended.  The  bivariate  Burr 
distribution  developed  in  this  paper  is  found  to  be  more  useful  in  application  than 
that  developed  by  Takahasi  (1965) . 
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